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Greek alphabet 


a A alpha b+. 466a p= a0 
BB beta kK K_ kappa aos + sigma 
y [ gamma A <A lambda 7. t= tau 

6 A delta uw M mu v  Y~ upsilon 
E E_ epsilon ea @ @® phi 

'- Z zeta ae eo A 

n H_ eta o QO - omicron wy W psi 

dv 0 theta f= pi w 2 omega 


2 SI units (Systeme International d’Unités) 


The SI system is based on six base-units, of which the following four are used in this 
course: 


length: the metre (abbreviated to m); 
mass: the kilogram (abbreviated to kg); 
time: the second (abbreviated to s); 


temperature: the kelvin (abbreviated to K). 


Of these, only the kelvin may be unfamiliar. To express a temperature in kelvins, add 
273.2 to temperature expressed in degrees Celsius or ‘centigrade’ (so the temperature of 
melting ice is 273.2K and that of boiling water is 373.2 K). The absolute zero of 
temperature is 0K. 


Units for other quantities are given in terms of the base-units, as in the following 
examples: 


velocity: metres per second (ms~*); 
acceleration: metres per second per second (ms~7”); 
angular frequency: radians per second (rads~'). 


Notice that ‘per’ is represented by an index —1, and that the usual laws of indices 


apply, so that ‘metres per second per second’ becomes ne Ss Ss 


The following combinations of base-units occur commonly enough to have special 
names: 


force: one newton (N) = _ one kilogram-metre 

per second per second (kgms 7); 
energy: one joule (J) = one newton-metre (kgm? s~?); 
power: one watt (W) |= = one joule per second (kg m?s~°) 


Other units are often expressed in terms of these subsidiary units, e.g. 


stiffness of a spring (see Unit 7): newtons per metre (Nm~?); 
dashpot constant (see Unit 8): newtons per (metre per second) (Nm's). 


To avoid very large or very small numbers we also use multiple and fractional units, for 
example the kilometre (km) which is 1000 metres. The most important prefixes for 
forming these are given below. 


giga = 1000000000 = 10° (abbreviated to G) 
mega = 1000000 = 10° (abbreviated to M) 
kilo = 1000 = 10° (abbreviated to k) 
milli = 1/1900 = 10°° (abbreviated to m) 
micro = 1/1000000 = 10~° (abbreviated to pz) 
nano = 1/1000000000 = 10-® (abbreviated to n) 


So, for example, the pressure of the atmosphere, which is about 10° Nm~?, is more 
conveniently written 100kN m~?, i.e. 100 kilonewtons per square metre. 


3 Algebra 


3.1 Polynomials 


A polynomial of degree n is an expression of the form 
@.8° +o, 16° as + 


where n is a positive integer, x is a variable and ao, a1,...,@, are constants with 
a, #0. A linear polynomial has n = 1, a quadratic polynomial has n = 2 anda 
cubic polynomial has n = 3. 


The roots of a polynomial p(x) are the solutions of the equation p(x) = 0. The number 
of distinct roots is (at most) equal to the degree of the polynomial. If a polynomial 
p(x) has a root a, then it has a factor x — a. If it has a pair of complex conjugate roots 
a and @, then it has a real quadratic factor (x — a)(x — @) = x? + bx +c, where 

b= —a—@ and c= aa. 


To factorize a polynomial is to express it as a product of two or more polynomials of 


lower degree. The factors of the difference of two squares x” — a” are (x — a)(x +); 


those of the perfect square x? + 2ar + a? are (x + a)(x +a), ie. (x +.a)?. 


Completing the square means writing a quadratic polynomial as a sum or difference 
of perfect squares, as follows: 


ax? + be +c=a(2+b/2a)* +k’, 


where k = ,/+(c — b?/4a), the sign being chosen to make k real. 


3.2 Summation notation 


If a1, @2,...,@n is a sequence of numbers, then 
Tm 
» ai means 0; + @9 + «*>+ @,=1 + Gp: 
i=1 
Further, if m is an integer less than n, then 


Tr 
ya Means. —-0,, 7:Gy,41 + °°? FF Oe ee 


ii 


3.3 Factorials and the binomial theorem 


For any positive integer n, we define n factorial, written n!, by 
ee (1) * 20. 
The first few are 1! = 1, 2! = 2, 3! = 6, 4! = 24. We also define 0! = 1. 


The binomial theorem states that, for any positive integer n, 


—] 
(a +b)" =a" +na”™ *b+ me nie +--+: 


! 
fee ee ee 
ri(n—r)! 


= n! | 
— ) ee ee 
er as 
In particular, 
(a+b)? = a* 4+ 2ad + O°, 


(a+b)? = a Sab 4 Sab + 8°, 
(a + b)* = a* + 4a°b + 6a7b? + 4ab® + 5. 
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3.4 Powers and logarithms 


If a > 0, the power function a® has the properties 


aera aa 
qgt/¥ = (at/¥)* is (a® 7", 
Ifa>0,a4l1andz > 0, the logarithm function log, x has the properties 
log, 1 = 0, 
log,(2~") = —log, , 
log, (xy) = log, x + log, y, 
log, (2/y) = log, « — log, y, 
log,(z") = y log, z. 
The power and logarithm functions are inverse functions. Hence 
¢ = log. y means <-y= a” 


and ia. means 4 = log, =. 


3.5 Partial fractions 

This technique is useful for breaking up functions of the form P(x)/Q(x), where P and 
Q are both polynomials, into a sum of parts which can be integrated easily. 

A simple case (see Unit 2) is the function 


ee ey 
a> ax? + ba+c’ 


where b? — 4ac > 0. The polynomial ax? + bx + c has two real roots: call them x; and 
XL. Then we have 


f(z) = 
We assume that numbers N; and N2 exist such that 


+ N N. 
a(z — 27 )(2—22) 2 22 


so that az + B= aN, (x — r2) + aNo(x — 71). 


axr+ B 
a(x — 21)(z — 22) 


Equating coefficients of x and the constant terms on both sides of the equation gives 


a= aN, a aNo, 
B= -azgl; — azz, 


a pair of simultaneous equations which can be solved for N; and N2. So we have 


f(z) = E+ 


e—ey- ££ oy 
which can be integrated using the table of standard integrals. 


The general procedure for partial fractions is as follows. 


$ 


Let the function to be simplified be f(x) = P(x)/Q(x), where P and Q are 
polynomials with real coefficients. 


If P(x) has degree equal to or greater than that of Q(x), bring it to the form 
P(x) = Pi(x)Q(x) + R(z), 
where P\(zx) is a polynomial and R(x) is a polynomial of degree lower than that of 
Q(z). Then 
R(x) 
Q(x) 


Factorize Q(x) into real linear or quadratic factors. 


f(x) = Pi(x) + 


For each non-repeated linear factor (x — x;) in Q(x), write a term N;/(x — x;) in 
the partial fraction expansion. 


For each non-repeated quadratic factor x? + b;z + c¢; (corresponding to a pair of 


————— mn the partial fraction 
x? +b;7+4+¢; 


complex conjugate roots in Q(a)), write a term 
expansion. 
For each repeated linear factor (x — x;)”* in Q(x), write terms 

(c—a2;) (x—2;) (2 — 2,)™ 
For each repeated quadratic factor (x? + b;x + c;)™, write terms 

(x? + bjx + c;) (x? + br + ¢;)™- 
Equate the sum of partial fractions to R(z)/Q(xz) and multiply both sides of the 
equation by Q(x). Then equate the coefficients of like powers of x on both sides, to 


obtain a set of simultaneous equations for all the unknown constants 
Dy, NG ee ae. ss 


3.6 Rules of algebra 


In the following summary of the rules or axioms of algebra, the symbols a, b, c stand 
for arbitrary elements of some set S (which can be the set of, for example, real numbers, 
complex numbers, square matrices, n x m matrices, or vectors) and the symbol * stands 
for the binary operation being considered, which may be addition, multiplication, etc. 


Definitions 

Closure: If a and b are elements of S, then so is a * b. 

Associativity: (a *b) #c = Ge ieee. 

Identity: There is an ‘identity’ element e such that axe =e*xa=a. 
(For addition it is called 0; for multiplication, 1 or J.) 

Inverse: To each element a there is a unique inverse element a 


such that a*#a = a%+a =e: 
(For addition @ is called —a; for multiplication, a~'.) 


Commutativity: axb=b*a. 


Distributivity: ax (boc) =(ax*b)o(axc). 


(over binary (For example, for addition and multiplication of real numbers 
operation o) a(b+c) =ab+ac.) 


8 


Summary of rules 


Multiplication 


Addition Real Complex | Square | Non-square | Vectors: | Vectors: 
(all cases) | numbers | numbers | matrices | matrices dot cross 
product | product 


. Closure 
. Associativity 
Identity 

. Inverse 

. Commutativity 
. Distributivity 
(over addition) 


Orr WN 


Notes 
(i) Every element except 0 has an inverse. 


(ii) All non-singular square matrices (i.e. those whose rows are linearly independent) 
have inverses. 


(iii) The product of an n x m and an r x s matrix (in that order) exists if and only if 
m =r, in which case the product is an n x s matrix. 


4 ‘Trigonometry 


4.1 Trigonometric functions and their inverses 


Definition for acute 
angles in terms of 
triangle above 


Range of Definition Range of values* 
values for x for inverse function 


*Sometimes called the principal value range of that inverse function. 


4.2 ‘Two useful triangles 


oe ES 


1 
From these it can be seen that, for example, sin : == 5 and tan > = /3. 


4.3 Trigonometric identities 


Pythagorean relations 
sin? 6+ cos? 6 = 1 

1 + tan? 6 = sec? 6 

1 + cot? 6 = cosec? 6 


Addition formulas 


sin(a + 3) = sinacos 3 + cosasin 3 
sin(a — 3) = sinacos GB — cosasin 3 
cos(a + 3) = cosacos 3 — sinasin 3 
cos(a — 3) = cosacos Z + sinasin 3 


tana + tan 
t = Sees 
ae) 1 —tanatanf 

t = 
aa ana — tan @ 


1+tanatan 
5 sin(a + 8) + $sin(a — 8) 
5 sin(a + 3) — §sin(a — B) 
cos a cos 3 = 4 cos(a + 3) + 5 cos(a — f) 
1 
2 


cos(a — 3) — $cos(a + 8) 


sin acos 3 = 


cos asin 3 = 


sina sin 3 = 


10 


Double angle formulas 


sin 2a = 2sinacosa@ 
cos 2a = cos’ a — sin 
2tana 


2q =1—2sin* a = 2cos?a—1 


ian 26 = > 
1 — tan’ a 
(1 — cos 2a) 


i 
2 
cos? a = $(1 + cos 2a) 


4.4 Polar coordinates 


The point P whose polar coordinates are [r, 6] has Cartesian coordinates (z, y), where 
g=rceos@, .y=7sing. 


The value of r is always positive (except at the origin, where it is zero). For a given 
point P the value of @ is not unique: we can add or subtract any integer multiple of 27 
and obtain another value for 6 which describes the same point. The value of 
satisfying —7 < 0 < 7 is called the principal value of 6, and is given by 


2 
arccos | ——— if y= a, 
fx? + y? 
§= 
Lt ; 
—arccos | —————— ee 
x? + y? 


5 Functions and graphs 


5.1 Functions 


Any expression or formula which involves a variable (say +), and whose value is 
uniquely determined by the value of z, is called a function of z. 


If a variable y is a function of z (i.e. is equal to a function of x) then we call x the 
independent variable and y the dependent variable, and we may write y = y(z). 
Here y(x) stands for the function of z (i.e. for the formula involving 2). 


A constant is a variable whose value remains the same whatever value we give to the 
independent variable. 


It is also possible to think of a function in a more abstract way, as a mathematical 
object consisting of three parts: 

(i) the domain, which is the set from which the value of z is taken, 

(ii) the codomain, which is a set in which the value of y will lie; 


(iii) the rule which associates to each x in the domain a unique element in the 
codomain, called the image of z. 
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If no domain and codomain are specified, then the domain is taken to be the largest 
possible set of real numbers for which the rule makes sense, and the codomain to be the 
set of all real numbers. 


If f and g are two functions, then their sum is a function f + g defined by 


(f +9)(x) = f(x) + (x) (for all x). 
Moreover, if A and B are any two numbers, then the function Af + Bg is defined by 


(Af + Bg)(x) = Af(x)+Bg(x) (for all z). 


The function f(g(x)) = f ° g(z) is called the composite function or composition for 
the functions f and g. 


The graph of a function f is the curve in the (z, y)-plane whose equation is y = f(z). 


A constant function is one that assigns the same image to all the elements in its 
domain. Its graph is a straight line parallel to the z-axis. A special case is the zero 
function, which assigns the image zero to all the elements in its domain. 


A linear function is one having the form az + 6, where a and 6 are constants. The 
constant a is the slope and the constant b is the y-intercept of the graph. 


A continuous function is one such that a small change in x cannot produce a large 

change in f(z) (i.e. the graph of f has no breaks or jumps). 

The limit at a of a function f, written lim f(x), is a number L such that we can make 
ra 


|f(a) — L| as small as we please by choosing any z sufficiently close to a. For a 
continuous function, lim f(z) = f(a). 
r—a 


The limit for large x of a function f, written lim f(z), is a number L such that we 
L— OO 


can make | f(z) — L| as small as we please by choosing any z sufficiently large. 


5.2 Graphs of some common functions 


See sketches on page 12. 


5.3 Curve sketching 


Suppose we wish to sketch the graph y = f(x), where f(x) is some given function. A 
possible procedure for doing this is given below. 


1. Check if f(x) is any standard function whose graph you already know or a simple 
combination of such functions. 


2. Determine how y behaves when z is very large and positive and when z is very 
large and negative. 


3. Find where the curve crosses the z- and y-axes, if at all. 


4. Look for any values of x at which f(z) is undefined, and examine the behaviour of 
f(x) near these values of z. 


5. Find if there are any local maxima, local minima or points of inflection in the given 
domain. Find these values of f(x) and the associated values of x. 


6. ‘Transfer the information found in Steps 3 and 5 to a sketch graph, and use what 
you know about the general shape of the graph to help you join the points with a 
smooth curve, taking into account any information found in Steps 2 and 4. If you 
are still unsure about any parts of the curve, choose suitable values of x and plot 
the corresponding points (x, f(2)) before completing the curve. 
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Sketch graphs of some common functions 
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6 Differentiation 


(Only ordinary differentiation is considered here. For partial differentiation, see the 
summary of Unit 25.) 


6.1 Notation and terminology 


If f is a function, its derived function or derivative f’ is defined by 


_ f(a +h) — f(z) 

/ —— 

a h | 

The process of calculating f’(x) from f(z) is called differentiation of f(x) with 


respect to x. Differentiation with respect to x is denoted by the symbol — written to 


the left of the expression or variable which is differentiated, so that, for example, 


d, 
A means f(z), 
d 
and if y is a variable that stands for f(x), then “Y is a variable that stands for [ge ee 


dx 


d 
To save space, we often print dy/dz in place of a 
x 


d 
When the independent variable is the letter 2, we often use a prime in place of —, so 


dx 
d ; 
that y’ means - When the independent variable is t (time), we often use a dot in 
lace of that ¢ = 
ace of —, so that WU means —. 
dt dt 


The derivative of a derivative is called a second derivative. For example, the second 
derivative of the function f, denoted by f”, is the (first) derivative of f’, defined by 


fiw +h) = f(a) 


7 oa 
2 
d*y d (dy 
tt ie ee : UNSES Pans 
The symbol y” or 728 used for = ( z), 
d? df -d 
and =e Se is used for = (F). 


Third and higher derivatives are defined and written analogously. The nth derivative 
of f is denoted by f(™). 


6.2 Rules of differentiation 


Sum rule: if u and v are functions of x, then 


d du dv 
aa +v)= as + — 
Constant multiplier rule: if k is a constant and u is a function of z, then 
d du 
Product and quotient rules: if u and v are functions of z, then 
d du dv 
ay hue) = roe 
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Chain rule or ‘function of a function’ rule: if f and g are two functions, then the 
derivative of their composition is given by the function of a function rule: 


d 

7, lf (o(@))] = F'(9(@))9°(2). 
A more common way to write this is the chain rule: 

dy _dydu 

dx dudz’ 
where (in this case) u stands for g(x) and y for f(u). 
Implicit differentiation: given an equation connecting two variables, say x and y, 
we calculate dy/dzx by differentiating both sides with respect to x and then solving the 
resulting equation algebraically for dy/dx (instead of solving for y before 
differentiating). 


6.3 Standard derivatives 


In each of the following cases the domain is the largest set of real numbers x for which 
the function is defined. 


Function Derivative 


z ax*-! (qa any number) 
log, x 1/x 

log.(—z) 1/2 

e* e* 

sin x COS x 

COS © —sinzx 

tan x sec? x 

cot x — cosec? xr 

sec x sec z tan z 

COSEC & — cosec © cot £ 


The following derivatives are given in a form that is useful for integration. The 
constant a is assumed positive, and in the derivatives of arcsec and arccosec it is 
assumed that x > a (rather than x < —a). 


Function Derivative 


1 


ee 


818 
Saisie 


arcsin ( 
arccos 


arctan 


Dae Re eee 


SiS gIi&®& eatxXg SI8 


arccot 


Q 

bo 
+;9 
8 

bo 


arcsec 


Pi CMM az nk 


la 
8 
8 

ss | 
| 
i) 

i) 


arccosec ( 


Qls 
Pediat 
8 
8 
wl] | 
l i) 
0) 
i) 


6.4 ‘Taylor series 


About zx = 0: 
ee 
t= * Se eo : 
oo ee 
eg ei 
sy = 
= oe a a 
2 3 
log.(1+2)=2- = 4 = =: (=o < 1), 
r(r —1) 


(a+2)"=a"+ra” *r+ 


e a’ *g*+--. (-a<a2<a) 


(where r is any real number and a is any positive number). 


In general, about x = 0 


f(x) = f(0)+2f'(0) + 


and about x =a 
f(a) = f(a) + (wa) f(a) + S(e—a)?f"(a) +>. 


The nth-order Taylor approximation to f(z) about x = a is 


i 1 
a Oty PO 


f(x) = f(a) +(e — a) f(a) + Hl — a)? f"(a) +--+ (ea) F(a. 
In particular, n = 1 gives the tangent approximation 

f(x) ~ f(a) + (x — a) f(a) 
and n = 2 gives the quadratic Taylor approximation 

f(x) ~ f(a) + (x — a) f'(a) + 3(x — a)" f"(a). 


These approximations are good when z is close to a. 
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7 Integration 


7.1 Standard integrals 


The following integrals are obtained by reversing the entries in the tables of standard 
derivatives. In each case n stands for a non-negative integer, a for a non-zero real 
constant and a for a real constant. 


Function Integral Domain 
grtt 
- +C All x 
n+1 
att 
r* (a#-1) 0 ae 
= (and x < 0 as well if a is an integer) 
1 
= es oe Se | 
: 
: log.(—r) +C xz <0 
1 
ee —¢@" £C¢ All x 
a 
1 
sin ax ——cosaxr+C All x 
a 
i 
cos ax —~sinar+C All x 
a 
1 
sec? ax —tanaxr +C a 
a 2 2 
1 
cosec? ax ——cotazr+C O<ar<T 
a 
1 T T 
sec az tanazr —~secar+C —— <ar< — 
a 2 2 
1 
cosecaxcotax ——cosecar+C O<ar<T 
a 
1 ae C 
Ss ae arcsin = + —< F< 
1 + 
: = — arctan = at; A 
a* +2 a a 
1 1 
2 gee (ee ae 
rv x2 — a2 a a 
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The following integrals may also come in useful. As before, the constant a must be 
non-zero, but b, a, 3,-y can take any real values (except where stated otherwise). 


Function Integral Domain 
1 1 
5 ~ loge(ax + 6) +C ax+b>0 
: Pies b)+C +b<0 
— log.(—ax — 
ax +b < a 
log, ax x(log.(ax) —1)+C ax > 0 
1 1 7 
t ——] C —— <a < — 
an ax Z og.(cos ax) + RES 5 
1 
cot ax — log. (sinaz) +C G<i0r< 7 
s : 
1 us 
sec ax — log.(secaxr + tanazx) + C 3 
a 
1 
cosec ax — log. (cosec ax — cotaxr) +C 2 ar< 
a 
1 1 at+zr 
== sq ee( $=) +0 In| < Jal 
1 1 t= 0 
—] C 
are 5, ee( ==") + > [al 
1 = 
5a lobe ( ste )+e x < —|a| 
2a =f 
1 
———— log.(@ + Vz? +a7)+C All x 
Ss g(r+V ) 
1 
aT log.(a + Vx? —a?)+C a > fa] 
x2 —a 
log.(—x — Vz? —a?)+C x < —|al 
1 1 
ee oe. EF + aa? + 25 +7) +C If ax? + 2624+ 7 
Van? + 28x +¥ va va has real roots, 
(a > 0) x > the larger root; if 
not, x is unrestricted. 
1 1 ar — G3 


x between the roots 
of —az?4+26r+7 


————— — —= arccos ———————. + C 
\/—-an? + 207% +7 Va BA +ay 
(a > 0,6? +ay>0) 
2 


Vz? +a? = log.(a + Va? +a?) + = Va? a? +C Allz 


Zz 
2 
xr? — a? ee ee xz > |a| 
2 
aomee + Vx? —a*)+ ae —-a*+C 2«<-—|a| 


2 


yar 2 pein = + Ve =? 40 |z| < |a| 
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For the following integrals, the domain comprises all real values of z. 


Function Integral 
eet 
e*" sin Ox Se (asin Gx — Bcos Br) + C 
+ 8 
Ped 
e** cas Oe ——,, (a.cos Br + Gsin Bx) + C 
az+ 3 
i 
xsinax ——cosazx + — sinar+C 
a a 
| ae : 
COS ax — sinaz + — cosar+C 
a a 
. x? 2c. 
xr sin ax —_— + — }cosazr+ —sinar+C 
a a? a? 
9 x 22 
r~ COS ax A sinax + — cosar + C 
a a a 
Tm 
: 4b nr = 
x2” sinax —=—=—(osdr + —~1 a” loosardr+C 
a a 
; i) Sees 
x2” cos ax ~~ sinaz —— | 2” 'sinardr+C 
a a 
: 1 : 
sin? ax —~¢ — —sin2ar+C 
2 4a 
1 i 
cos’ ax —7 + —— gen ears G. 
pe 4a 


1 
sinaxcosaxr — cos 2axr+C 
a 


In the following three cases we require a £ 3 and a F# —{. 


Function Integral 
sin(a—fB)x  sin(a+ B)z 
sinazsin Be =£————— — ———_ + C 
me “aa—p) a+) 
sin(a—)x_  sin(a+ f)zx 


cos ax cos Bx 


%a—p) | at B) 


sinaxcosfx — Sia wars ic = = < a sé 


2(a — f) 


7.2 Rules of integration 


1. The fundamental theorem of calculus 


[oar = F(x)+C 
dF (x) 
a 


if and only if f(z) = 


2. Linearity rule 


[tos(e) + 69(2) de fc 


=a f(e)ar +8 fg (x) dz. 
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3. Substitution rules 


For integrals of the form To find / sin(ax + b) dx = Ih, let 
r= foe aes u = az + b; then adz = du. So 
- | d I : / in(az + b)ad 
= — | sin(az a dx 
substitute u = f(x) and use - dz = du to — 


1 1 
obtain ~ | sinudu = ——cosu+C 
= —-—cos(az+ b)+C. 

a 


Otherwise, find a suitable relation between u To find 2 == J5, let 
(1 + x?)3/2 


and x and use 
x = tanu; then dz = sec” udu. So 


dx 
de = 2, du. = / sec? u du 
eS (1 + tan? u)3/2 


sec? udu 
sec? u 
a [cosudu =sinu+C 


= sin(arctan x) + C. 
4. Integration by parts 


Let u = u(x) and v = v(z); then To find [ cosa ee ee 
[2 di ages [eS dz. and v = sinz; then 
dx dx d 
i= [eZ bing dx 
= ais = [sine x ldzx 


= ¢sinz + cosz+C. 


7.3 Integration methods 


Rational functions 

To integrate a rational function, that is, one of the form P(x)/Q(x) where P(x) and 
Q(x) are polynomials, we first express the function as a sum of partial fractions (see 
Subsection 3.5). Then integrate each partial fraction as follows: 


d 
/ = 7 a standard form; 
ax 
d 
/ ae substitute u = ax + b; 
ax 
/ a dx with b? < 4ac, 
(ax* + be +e) substitute 2az + b = V/4ac — b? tanu. 


Irrational functions 


Integrals containing the expressions Va? — x7, Vx? — a”, Vx? + a? can often be 


simplified by the following substitutions: 


for (a2 —2* use x=asinu; 
for Vz —a* use x =asecu; 
for J 2? + a2 use x=atanu. 

More generally, for \/ax? + Gx +7: 
ifa>Oand @? <4ay use 2ar+GB=/4ay— @ tanu:; 
ifa>Oand @?>4ay use 2ar+B= V8 —4ay secu; 
ifa<Oand B*>4ay use 2ar+3= VV? —4aysinu. 


(In the remaining case the square root is always imaginary.) 
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Exponential and trigonometric functions 
gee 
For ] P(x) ¢ sinax } dz, where P(z) is a polynomial, integrate by parts (several times 
Cos ax : 
Vv 


if necessary) with u = P(z), = e*” or sinaz or cosaz. See also the ‘standard 
xr 


integrals’ table for / x” sinax dx and J x” cosax dx. 


For / pees = eS : dx, see the ‘standard integrals’ table. 
cos Bx 


For / sin” az cos" ax dz, 


ifnmisodd use u=sinaz; 
ifm isodd use u=cosaz; 


if m,n are both even, then 
oi 0 se 0 = faa, 
im <( wee w= cote7, 
ifm,n both >0 use sin?az = $(1—cos2az), cos*az = $(1+ cos 2az). 


For other integrals involving trigonometric functions only, it may be helpful to express 
them in terms of sinaz and cosaz and then use the substitution t = tan an, according 
to which 

2t 1—¢?? 2 dt 


T4+f’ ee a ee car= Soaeys 


anar= 


ee a 
2t 


5 


For integrals involving e*” only, the substitution u = e®” may be useful. 


§ Unit outlines 


Unit 1 Recurrence relations 


Section 1: Introduction to recurrence relations 


1. An equation relating each term in a sequence to one 
or more previous terms, which can be used to generate 
the whole sequence once the first few terms are known, is 
called a recurrence relation. 


2. A recurrence relation is of kth order if the difference 
between the highest and lowest subscripts in it is k. To 
generate a sequence using such a recurrence relation we 
need in addition the values of the first k members of the 
sequence; any set of conditions from which they can be 
deduced are called the initial conditions. 


3. A recurrence relation is linear if it can be written as 
Ur+1 = ArUr + bp Ur—1 + Crlp—2 1° ** -F De, 
where p, and the coefficients a, b;, cr, ... may depend 


on r, but do not depend on any of the us. Otherwise it is 
said to be non-linear. 


4. A linear recurrence relation has constant 
coefficients if none of the coefficients a,, b;, cr, ... 
depend on r. (But p, may depend on r.) 


5. A linear recurrence relation is homogeneous if 
pr = 0 for all values of r. Otherwise it is said to be 
non-homogeneous (or inhomogeneous). 


6. The general solution of a recurrence relation is a 
formula containing one or more arbitrary constants, from 
which all the particular solutions can be obtained by 
giving the arbitrary constants particular values which 
may be determined from the initial conditions. 


7. The general solution of the linear 
constant-coefficient first-order recurrence relation 


Ur+1 = aUr 7 P 
is 
ee a-—1l 
A+ np 


where A and B are arbitrary constants. 


Section 2: Linear second-order recurrence 

relations 

To obtain the general solution of the linear homogeneous 

constant-coefhicient second-order recurrence relation 
Ur+1 = aUr + bury_1, 


we first consider the auxiliary equation 


xz? =ax +b. 


There are three cases to consider. 


(i) Ifa? + 4b > 0, the auxiliary equation has two 
distinct real solutions, A and yp, and the general 
solution for the recurrence relation is 


tin = AX + By. 
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(ii) If a? + 4b = 0, the auxiliary equation has a repeated 
real solution 4, and the general solution is 


Un = (A+ Bn)X". 


(iii) If a* + 4b < 0, the auxiliary equation has no real 
solutions. (See Unit 5, Section 5.) 


Section 3: Numerical difficulties 


1. If Z is an approximation to some quantity whose true 
value is xz, we call  — x the absolute error in z, and we 
call (¥ — x)/zx the relative error in z. 


2. Consider any problem in which we are given numerical 
data and we want to calculate some result from them. 
The problem is said to be absolutely ill-conditioned 
(with respect to small changes in the data) if a small 
absolute error in the data gives rise to a significantly 
larger absolute error in the result. It is said to be 
relatively ill-conditioned if a small relative error in 
the data gives rise to a significantly larger relative error 
in the result. It is said to be absolutely 
well-conditioned if the absolute error in the result is 
smaller than the absolute error in the data, and 
relatively well-conditioned if the relative error in the 
result is smaller than the relative error in the data. 


3. The scale factor of a calculation is the factor by 
which a small error in the data is multiplied to give the 
corresponding error in the result. If the magnitude of the 
scale factor is significantly greater than 1, then the 
problem is absolutely ill-conditioned with respect to small 
changes in the data; if the magnitude of the scale factor 
is less than 1, then the problem is absolutely 
well-conditioned. 


4. Consider the problem of calculating u, from uo, using 
the recurrence relation 

Ur+1 = Qty +p. 
The scale factor for this problem is a”. If the magnitude 
of a” is less than 1, the problem is both absolutely and 
relatively well-conditioned with respect to small changes 
in uo. If the magnitude of a” is significantly larger than 
1, then the problem is absolutely ill-conditioned. For any 
a such that |a| > 1, the problem is therefore absolutely 
ill-conditioned for sufficiently large n. It is relatively 
ill-conditioned if, in addition, the condition 


is satisfied. 


5. If |a,| > 1 for all r, then the problem of calculating un 
from uo using the more general recurrence relation 


Ur+1 = ArUr + Pr 
is absolutely ill-conditioned (with respect to small 
changes in uo) for sufficiently large n. It is relatively 
ill-conditioned if, in addition, the initial value uo is close 
to a value for which u,, the true solution of the 
recurrence relation, remains small for large values of n. 
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Section 4: Living with ill-conditioning 


1. Mortgage problems, relating to a situation where £X 
are borrowed at 1% per annum interest and repaid by 
monthly payments of £M, can be modelled using the 
recurrence relation 


Urg1 = (14+ 1/100)u, — 12M, 


where u, is the amount owed at the end of the rth year. 
The initial condition is uo = X. 


Mortgage problems are often highly sensitive to small 
changes in the data, particularly the values of J and M. 


2. Sometimes ill-conditioning can be cured by a 
reformulation of the problem. For example, in the 
recurrence relation 

Ur+1 = GrUr + Pr; 
with |a,| > 1 for all values of r, the solution which 
remains small cannot be reliably calculated by forward 
recurrence, i.e. starting with wo and calculating wi, 
u2, ... in succession, because this problem is relatively 
ill-conditioned. The reformulation which cures the 
ill-conditioning is to use backward recurrence instead. 
This means rewriting the recurrence relation with uy as 
the subject of the formula, 


Ur+1— Pr 
up = ———— 
ar 
and calculating values of u, from it working down from a 
large value of r, say R. Even though we know very little 
about the value of ur, it is easy to obtain accurate 


results because the new problem is very well-conditioned. 


Unit 2 Differential equations I 
Section 1: Direction fields 


1. A first-order differential equation is an equation 
connecting the derivative dy/dz of some unknown 
function to the variables x and y which are related by 
that function. By solving the equation algebraically for 
dy/dz we can bring it to the standard form 


oH = m(2,9), (1) 
where m(z, y) stands for some formula which may 
contain either or both of the variables x and y. 


2. By a solution of a differential equation we mean a 
function—that is, a rule or formula giving y in terms 
of z—which satisfies the differential equation. 


3. A direction field is a rule or formula which specifies 
a slope at each point in the (z, y)-plane (or some part of 
it). To each direction field there corresponds a differential 
equation, constructed by taking m(z, y) in (1) to be the 
formula giving the slope of the direction field at the point 
(x,y). If y = f(z) is the equation of a trajectory of the 
direction field—that is, a smooth curve at every point of 
which the slope of the curve matches that of the direction 
field—then y = f(z) is also a solution of the 
corresponding differential equation. The direction field 
has an infinite family of trajectories, each one 
corresponding to a different solution of the differential 
equation. 


4. For some differential equations, an exact formula can 
be found for the solutions. This formula, called the 
general solution, contains an arbitrary constant of 
integration, often denoted by C. Different numerical 
values for C give different solutions to the differential 
equation (usually called, in this context, particular 
solutions to distinguish them from the general solution) 
and therefore label different trajectories of the direction 
field. 


5. Whether or not it is possible to find a formula for the 
solution of a differential equation, the direction field is a 
useful source of information about its solution; in 
particular it can be used to sketch solution curves 
(trajectories) and so obtain an idea of their qualitative 
behaviour. 


Section 2: A numerical method 


Given a differential equation 


aye 
ae = m(x, y) 


and the condition that y(xzo) = yo, where xo and yo are 
specified, then Euler’s method can be used to find y(6), 
where 6 is a given value of z, as follows. 


(i) Choose h = (b— zo)/n, where n is the number of 
steps to be used. 


iy otra 
Y,41 = ¥o+ bones ¥,), 


where xz, = 20 + rh, so that rz, = 6, and Yo = yo. 


— 1, use the recurrence relation 


(iii) Then Y,, is an approximation to y(b). 


The error in the approximation to y(b) is roughly 
proportional to h, and therefore inversely proportional to 
the cost (i.e. number of steps) of the calculation. 


Sections 3 and 4: Exact methods 


1. To look for an exact solution to a given differential 
equation, first bring it to the standard form 


ot = m(z,y). 
Then 
(i) if m(z,y) has the form 
m(x,y) = f(x)g(y) 


(where either or both of f(x) and g(y) may be 
constants) use separation of variables (see below); 


(ii) if m(z,y) has the linear form 
m(x,y) = k(x) + U(a)y 
(again either or both of k(x) and I(x) may be 


constants) use the integrating factor method (see 


below). 


2. Separation of variables (Procedure 3.3) 


To solve the differential equation 


d 
7. =F (z)9(y) (gly) #9): 
= fe sides by g(y) and integrate, to obtain 


le dy = f $2) )dx + C. 


Perform - two integrations and make y the subject of 
the resulting equation. 


3. Integrating factor method (Procedure 4.2) 
To solve the linear equation 


oY = ko) + Un)y, 
subtract l(x)y from both sides and then multiply both 


sides by the integrating factor 


a | / (-U(2)) de | 


The resulting equation has the form 
dy | dp(x) _ 
pet eee) 
and integrates directly to give 


sey = / ieee 


which is easily solved for y to give the general solution. 


4. Given the general solution of a differential equation 
involving an arbitrary constant, say C’, we can determine 
the particular solution satisfying a given further 
condition by substituting the general solution into the 
further condition. This gives an algebraic equation which 
we can solve for C’. Inserting this value of C into the 
general solution of the differential equation gives the 
required particular solution. 


5. If you need help with the integrations, consult the 
Integration section of this Handbook. 


Unit 3. Animal populations: their 
growth and exploitation 


Section 1: How to model populations 


Mathematical modelling is the use of mathematics in 
solving real-world problems. Empiricism means basing 
one’s model on experiment, without any theory; in the 
study of population changes it means fitting a simple 
mathematical relation to the data without seeking to 
explain why it should have that particular form. If the 
empirical relation is then applied outside the range 
covered by the original data, we are extrapolating, a 
notoriously inaccurate procedure. 


A useful tool for setting up models which do have some 
theoretical basis is the input—output principle: 
increase in population (i.e. accumulation) 
= input — output. 
The exponential model is obtained by assuming that 
(i) the population P(t) can vary continuously, 
(ii) migration is negligible and (iii) the numbers of births 
and of deaths in any short time interval of duration h are 
respectively BP(t)h and CP(t)h, where B and C, the 
proportionate birth rate and proportionate death 
rate, are taken to be constants. The resulting equation is 
dP 
— (B-—C)P. (1) 
The long-term predictions of this model are 
(i) unbounded growth if B > C, (ii) decay to zero if 
B < C and (iii) constant population if B = C. Case (i) in 
particular is clearly unrealistic, but for some populations 
these predictions fit the real data reasonably well over 
limited time intervals. 
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Section 2: The logistic model 


In the logistic model we assume that B and C, instead 
of being constants, depend on the total population in 
such a way that B — C, the proportionate growth 
rate, decreases linearly with the population P. The 
resulting equation is the logistic equation 

ae ~*(1- 9)? 2 
where a and M are constants. Built into this model are 
the assumptions that 


(i) there is an equilibrium population level, M; 


(ii) the proportionate growth rate depends only on the 
instantaneous value of P; 


(iii) the proportionate yrowth rate decreases linearly with 
P, being given by a(1 — P/M). 


The solution of the logistic equation is 


M 
ly Sees 
(t) 1+ (M/P(0) — 1)e-2’ 
and the graph of this solution (for 0 < P(0) < M) is 
called a logistic curve. 


t 


Two ways of comparing this model with observed values 


for P(t) are: 


(i) estimate M from the property that P(t) is close to 
M for large t, then plot log.(M/P(t) — 1) against t; 


(ii) if data are available at equally-spaced instants 
fot, plot (Prat eee coeainast. Fy). 


In each case the logistic model predicts that the graph 
will be a straight line. 


Though far from being a ‘law of nature’, the logistic 
model agrees quite well with data for some populations, 
particularly under laboratory conditions; but for wild 
populations the agreement is more likely to be qualitative 
than quantitative. In cases where the model does not fit 
the data well, it is usually possible to see that one or 
other of the assumptions (i) to (iii) above is not realistic. 


Section 3: Fishing 


The following procedure, known as the Schaefer model, 
can be used to estimate optimum fishing effort and 
maximum sustainable yield for a fishery. 


(i) Plot yield per unit effort (a measure of population) 
for each year against the mean fishing effort over the 
preceding few years. 


(ii) Fit a straight line to the points plotted in (i) and 
assume this to give the relation between the 
long-term sustainable yield per unit effort and 
sustained fishing effort. 


(iii) Deduce the relation between sustainable yield and 
sustained fishing effort, and hence find the sustained 
fishing effort that gives the maximum sustainable 
yield. 
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The Schaefer model assumes that 
(i) the modified logistic model (3) below applies; 


(ii) the historical variations of fishing effort were slow 
enough for the fish population in each year to be 
close to the equilibrium level appropriate to the 
mean level of fishing effort exerted over the previous 
few years; 


(iii) catch is proportional to the product of fish 
population and fishing effort. 


Section 4: Models of exploitation 


We can modify the logistic model to include the effect of 
exploitation by human beings. The new equation is 
dP - 
—=a(1—-—]P-C 3 
a ( 3 | (3) 
where C is the rate at which members of the population 
are caught, i.e. it is the yield. Two exploitation 


strategies are considered. 


(i) Constant catch, i.e. C = constant. Here the 
conclusions are 


(a) if C exceeds the maximum sustainable yield 
aM/4, the population will dwindle to zero; 


(b) even if C < aM/4, the population will dwindle 
to zero if P(0) is too small; 


(c) if C <aM/4 and P(0) is large enough, then for 
long times P(t) will approach an equilibrium 
value smaller than M. 


(ii) Constant effort, i.e. C proportional to P, say 
C = \fP, where f (the effort) and 4 are constants. 
Here the conclusions are 


(a) if \f >a, then the population will dwindle to 
ZeTO; 


(b) if \f <a, then the population is viable whatever 
the value of P(0); 


(c) in the long run, maximum yield is obtained with 
a value of Af = 5a. 


A more general model, in which (3) is replaced by 
aP 
—=f(P)-C 


with f(P) a suitable function of P, gives qualitatively 
similar conclusions. The qualitative features of the above 
conclusions are therefore more reliable than the 
quantitative ones. 


An important difference between the two strategies is 
that if we aim for maximum yield, the constant catch 
strategy is unstable—a small variation in the conditions 
can cause a large change in the result, and can in this case 
cause the population to dwindle to zero—whereas the 
constant effort strategy is stable—a small variation in 
the conditions produces only a small change in the result. 


Unit 4 Newtonian mechanics in one 
dimension 


Section 1: The description of motion 


1. A particle is a material object whose size and 
internal structure are negligible, so that at any given time 
it is located at a single point. 


eS 


2. The motion of a particle along a straight line is 
described by a position function z(t). For each instant 
t in time, this function specifies the position x(t) 
occupied by the particle. 


3. The velocity v of the particle is the derivative of the 
position function with respect to time, that is, 
_ dz 
Se 


4. The acceleration a of the particle is the derivative of 
the velocity function with respect to time. In symbols, 


dv d*x au 


i a a ae 


5. In the special case of uniform motion, by which is 
meant motion with constant velocity vo, we have 


= wet vot; 
VU = V0; 
a0. 


where Zo is the position at time t = 0. 


6. In the special case of motion with constant 
acceleration ao, we have 


2=270+ vot + Laot”, 


v=v90+ aol,—7 


v* = vo + 2ao(z — 20), 


where zo and vo are the position and velocity of the 
particle at the instant t = 0. 


Section 2: Newton’s laws of motion 


1. Newton’s first law states that, in the absence of a 
force, a particle either stays permanently at rest or moves 
at constant velocity. 


2. The force acting on a particle at a given instant is 
represented by a single number, F’, called the force’s 
gz-component. The force has magnitude |F'|, and is in 
the direction of increasing z if F > 0 or in the direction 
of decreasing zx if F < 0. Often ‘the force with 
xz-component F” is shortened to ‘the force F”. 


3. Newton’s second law states that if a particle of 
mass m experiences a net force with z-component F’, 
then its acceleration a is given by Ff’ = ma. 


4. If more than one force acts on a particle, then the net 
force or total force F' to be used in Newton’s second 
law is the algebraic sum of the individual forces, that is, 
the sum of their x-components. 


Section 3: First model for the motion of a 
falling object 
In this section Newton’s second law is used to model the 


motion of an object falling under the force of gravity, 
where the effects of air resistance are ignored. 


1. Newton’s second law F’' = ma may in principle be 
integrated by putting 


Sa < (giving v in terms of t) if F = F(t) or 
F = F(v); 


(i) a= ve (giving v in terms of x) if F = F(z) or 


Fe =F), 


2. Any object near the Earth’s surface is pulled 
downwards by the force of gravity. The magnitude of 
this force is proportional to the object’s inertial mass m, 
that is, 


force of gravity = mg downwards, 


where the constant of proportionality g is approximately 
9.81ms ”. 


Section 4: Second model for the motion of a 
falling object 


In this section Newton’s second law is used to model the 
motion of an object falling under the forces of gravity and 
air resistance. 


1. The motion of any object through the Earth’s 
atmosphere is opposed by the force of air resistance, 
whose magnitude R depends on the object’s speed |v], 
shape and size. For small speeds it is appropriate to use 
the linear approximation 


R ~ ki|vl, 
whereas for larger speeds the quadratic approximation 
R ~ kolo? 
is appropriate. Here k; and kz are constants which 
depend on the object’s shape and size. 


2. The constant speed |vp| at which an object can fall is 
called the terminal speed. Using Newton’s first law, 
this occurs when R = mg. 


Unit 5 Complex numbers 


Section 1: Introducing complex numbers 


The natural numbers are 1, 2, 3, ...; the integers are 
these together with their negatives and zero; rational 
numbers are numbers of the form m/n, where m, n are 
integers (n 4 0); real numbers can be thought of as 
possible lengths of line segments and are often 
represented geometrically as points on a number line. 
An irrational number is one that is real but not 
rational. 


A complex number can be defined as a pair of real 
numbers (xz, y) obeying the following laws of addition and 
multiplication: 

(v1, y1) + (x2, y2) a (x1 =e ete y2), 

(21,41) X (@2, y2) = (w1@2 — yiy2, Tiy2 + yi@2). 
Like the rational and real number systems, the complex 
number system is arithmetically complete in the sense 
that the sum, difference, product and quotient of any two 


numbers in the system (with the single exception of 
division by zero) is again a number in the system. 


If z is a complex number (2, y), then the real part of z is 
Rez = z and the imaginary part is Im z = y. Two 
complex numbers are equal if and only if their real parts 
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are equal and their imaginary parts are equal. A complex 
number is (purely) real if its imaginary part is zero, and 
(purely) imaginary if its real part is zero. 


The following notation, called the Cartesian form, is 
customary for complex numbers: 


(x,y) is written x + iy, 

(x, 0) is written z, 

(0, y) is written zy or yi, 

(0,1) is written 2, 

(0,0) is written 0. 
Calculations with complex numbers are then easily 
carried out by following the usual rules of algebra with 


the additional rule i? = —1. A list of these rules is given 
in Subsection 3.6 of this Handbook. 


The complex conjugate of a complex number 
z=2-+1y is the complex number Z = x — iy, which has 
the property that zZ is real and equal to x” + y?. To 
calculate the quotient of two complex numbers, multiply 
numerator and denominator by the complex conjugate of 
the denominator. 


Section 2: Geometrical representation 


The Argand diagram is a representation of the complex 
number system, in which the complex number z + ty is 
represented by the point in the plane with Cartesian 
coordinates (xz, y). Addition of complex numbers can then 
be represented by a parallelogram construction, and 
complex conjugation corresponds to reflection in the 
X-axis. 


Section 3: Polar form 


If the point representing a complex number z has polar 
coordinates [r, 6] (see Subsection 4.4 of this Handbook), 
we can write z in the polar form 

z=r(cos@+isin@). 
The modulus of z is |z| = r, and 6 is called an 
argument of z. If @ lies in the principal value range 
—mt <@< 17, we call 6 the argument of z and write 
6 = Argz. In terms of zx and y, the real and imaginary 
parts of z, we have 


jz] = V2? +y?, 
arccos(z/|z|) if y > 0, 
Arg z = 
—arccos(z/|z|) if y <0. 
When two complex numbers are multiplied, the moduli 
are multiplied and their arguments are added. In 
particular, de Moivre’s Theorem states that, for any 
integer n, 


(cos@+ isin 6)” = cosné + isin né. 


Section 4: Euler’s formula 


Euler’s formula 

e’’ = cos + isin 
enables us to write a complex number z in exponential 
form as 

id 

re. 
where r = |z| and 6 = Arg z, and to define the 
exponential function of z as 


xr+iy 


exp z =e” =e = e"(cosy + isiny), 


where z = Rez and y = Imz. 
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Euler’s formula also gives the results 


1 i 4 ‘ 1 a —1 
cos8 = =(e" +e rs sind = = (e" —e ys 


which can be used to derive trigonometric identities. 
A sinusoidal function of time t is one having the form 
f(t) = Acos(wt + $), 


where the constant A is called the amplitude, w the 
angular frequency and ¢ the phase. The function 
repeats after a time T' = 27/w, called the period. The 
complex number Ae’ is called the phasor of the 
sinusoidal function. Sinusoidal functions with the same w 
can be added by adding the corresponding phasors. 


Section 5: Complex roots 


If a = Re’® is a complex number and n is a positive 
integer, the nth roots of a, written a'/” are the 
solutions of the equation z” = a. There are precisely n of 
them, given by 


ai/” = Ri" exp (2+) (k=0; 1. ae. 


The fundamental theorem of algebra states that any 
polynomial of degree n, 
P(z) = an2z” + an—12" > +++ +a1z+ 40 
(with an # 0), can be factorized into n factors of degree 1, 
An(z — 21)(z — 22)... (2% — 2n), 
where z1,...,2n are complex numbers. These numbers 
are called the roots of the polynomial and are the 
solutions of the equation P(z) = 0. If the coefficients 
Qn,Qn—1,---, 40 are all real, then the non-real roots (if 
any) occur in complex conjugate pairs. 


For second-order recurrence relations of the form 
Ur+2 = QUr+1 a bur, 


if the auxiliary equation z* = az + b has non-real 
solutions re*’’, the general solution is 


un = 7" (Acosné+ Bsinné). 


Unit 6 Differential equations I] 


Introductory definitions 


A second-order differential equation is an equation 
relating the second derivative of an unknown function to 
one or more of the following: known functions, the 
unknown function, and the first derivative of the 
unknown function. The general solution of such an 
equation normally involves two arbitrary constants. If a 
second-order differential equation can be written in the 
form 
2 

p(x) <4 + a(2) 4 + r(ay =f) (1) 
where p, g, r and f are given functions (and p(z) is not 
the zero function) then it is linear. If in addition f is the 
zero function, then the equation is homogeneous; 
otherwise it is inhomogeneous (or 
non-homogeneous). If a linear second-order differential 
equation has the form 

2 

at - pe ey = F(z), (2) 
where a, b, c are constants with a # 0, then it is said to 
have constant coefficients. 


Section 1: Homogeneous equations 


To solve a homogeneous second-order linear equation 
with constant coefficients, i.e. one of the form 

d*y | ,dy 

—= + b— = 3 
apply the following procedure. 


(i) Write down the auxiliary equation 
ad? + bX be 8, (4) 
and solve it algebraically for A. 


(ii) (a) If (4) has two real distinct solutions A; and 2, 
then the general solution of (3) is 


Woy = Ae*!* + Be”, 
where A and B are arbitrary constants. 


(b) If (4) has only one real repeated solution A; (i.e. 
ad? + b\ + c is a perfect square), then the 
general solution of (3) is 


y(z) = (A+ Ba)e!”. 


(c) If (4) has complex conjugate solutions 
A1 =a+iB and A\2 = a — i, then the general 
solution of (3) is 


y(x) = e**(Acos Gz + Bsin Gr). 


If the differential equation is homogeneous but does not 
necessarily have constant coefficients, the following 
theorems tell us something about the solutions. 


Theorem 1.1: If the functions y; and y2 are solutions of 
a homogeneous linear differential equation, then any 
linear combination of y; and yp is also a solution. 


Theorem 1.2: If, in addition, the differential equation is 
of second order, i.e. of the form 
2 

p(x) + a(x) +r(x)y = 0 
with the coefficient p(x) never taking the value zero, and 
if the solutions yi and y2 are linearly independent, then 
the general solution of the differential equation is 
Ay: + Byz2, where A and B are arbitrary constants. 


By a linear combination of two functions yi and y2 
with the same domain we mean a function of the form 
c1y1 + c2y2, where c; and c2 are constants. The functions 
y1 and y2 are said to be linearly independent if neither 
is a constant multiple of the other. 


Section 2: Inhomogeneous equations 


For the general inhomogeneous Equation (1), we define 
the associated homogeneous equation to be 


p(x) + q(x) + r(x)y = 0. (5) 


The general solution of (5) is called the complementary 
function of (1), and is denoted here by y-. The general 
solution of (1) is (Theorem 2.1) 


Y = Vc T Up, 
where y, is any particular solution of (1). 


For certain constant-coefficient inhomogeneous equations, 
particular solutions can usually be found by substituting 
the function for yp suggested below into the differential 
equation and choosing the constants m, n or z (complex) 
so that the equation is satisfied. 


constant 


mztn 
ar 


constant 
kz +l 


ke” me 
kcoswx+lsinwx | mcoswx+nsinw2z 
or Re(k —il)e’*” | Reze*** 


In the above, k, 1, a and w denote given constants. 


Exceptionally, the differential equation may have no 
solution of the suggested form; in that case, multiply the 
suggested function by z and try again. 


The principle of superposition: If u(x) is a solution of 
Equation (1), and if v(x) is a solution of the equation 
obtained by replacing f(z) in (1) by g(z), then 

au(x) + bv(xz) (with a and 6 constants) is a solution of the 
equation obtained by replacing f(x) in (1) by 

af (x) + bg(z). 


Section 3: Initial and boundary conditions 


In an initial condition problem for Equation (1), the 
values of y and its derivative are both specified for the 
same value of z. Such a problem always has a unique 
solution (so long as p(x) never takes the value zero). 


In a boundary condition problem, one condition is 
placed on the value of y (or its derivative) at each of two 
different values of x. Such a problem may have a unique 
solution, or no solution, or an infinite number of 
solutions. 


Section 4: Nature of solution functions 


For the homogeneous Equation (3), with a and c positive, 
the solution y(x) may depend on z in various ways, 
depending on the value of b: 


b < —2,/ac 


(for large x) increasing exponential; 


—2,/ac<b<0_ oscillation with increasing amplitude; 


b= 9 sinusoidal oscillation; 
0<b< 2,/ac 
2,/ac <b 


For the inhomogeneous equation 


2 
le + pty +cy = kcoswxr+lsinwaz, 
fe dx 

with a, b, c, k, 1, w constants and a, b, c all positive, the 
complementary function decays to zero as x increases, 
and is called a transient. For sufficiently large x, any 
solution of the equation is close to the unique purely 
sinusoidal particular solution, which is called the 


steady-state solution. 


(for large x) decreasing exponential. 


Section 5: A numerical method 


Any second-order differential equation 


se (« = 
dx2 ie 


can be written as two simultaneous first-order equations 


dy 
Se 
<= = m(z,y,2). 


Euler’s method approximates these by the recurrence 
relations 


oscillation with decreasing amplitude 
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Pega — Ys a hZ,, 
eae = Ze 4- hm(2,, ‘a bee 


where h is the step size, 7, = 29 + hr and Y,, Z, are the 
approximations to y(z,) and 2z(z;,). 


Unit 7 Oscillations and energy 


Section 1: A home-made oscillating system 
Section 2: The perfect spring 


In this section a perfect spring is used to model vibrating 
systems. 


1. A perfect spring is characterized by two positive 
constants: its natural length and its stiffness. It obeys 
Hooke’s Law, exerting the following force on any object 
attached to either of its ends: 


(i) when the spring is extended, the force is directed 
towards the centre of the spring and has magnitude 
(called the tension) equal to stiffness x extension; 


(ii) when the spring is compressed, the force is directed 
away from the centre of the spring and has 
magnitude (called the thrust) equal to stiffness x 
compression. 


2. Ifa particle in one-dimensional motion is attached to 
perfect springs then the equation of motion of the 
particle can be derived by considering just one 
configuration in which each spring is designated as being 
either extended or compressed. 


3. The equation of simple harmonic motion is 
f+w*r= 0, 

whose general solution can be written either as 
x = Bcosuwt+ Csinwt 

or as 
x = Acos(wt + ¢). 


The constant w is called the angular frequency and the 
constant ¢ the phase angle of the oscillations. The 
period of the oscillations (time for one complete cycle) is 
T= 27/w. The frequency (number of cycles per second) 
is f = 1/r=w/2n. The amplitude of the oscillations is 
A. The constants A, ¢ and B, C are related via the 
equations 

A= 7 B*+C? 
and 

= arccos(B/A) C <0, 
Re ee arccos(B/A) C > 0. 


4. For a particle in simple harmonic motion, the 
equation ¢ + w?x = 0 describes the motion when the 
origin of the z-axis is chosen to be at the equilibrium 
position of the particle. If another origin is chosen then 


the equation of motion becomes 


+s 2 2 
Lia 2S we, 


where x = Z- is the particle’s equilibrium position. For a 
particle of mass m attached to one end of a perfect spring 
of stiffness k whose other end is fixed, the angular 
frequency w is given by w* = k/m for either vertical or 
horizontal motion. Here gravity is taken into account in 
the vertical case, but other forces are ignored. 
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Section 3: The conservation of energy 


1. The relationship between a force with z-component 
F(x) and the potential energy function U(z) is 


iar - | Fede, 
or, equivalently, 


dU 


The point at which the potential energy function is zero 
is called the datum. This can be chosen to be any 
convenient point. 


2. In particular, the gravitational potential energy is 
U = mg x height, 


where the height is measured from the chosen datum. 


3. The potential energy of a perfect spring is 


U = 5 x stiffness x (extension). 
4. The law of conservation of mechanical energy 
states that if the total force acting on a particle depends 
only upon the particle’s position then the quantity 


imv’ +U(x) =E 
remains constant throughout the particle’s motion. The 


quantity 

‘a imo" 
is known as the kinetic energy of the particle, and the 
constant EF as its total mechanical energy. 


5. The region of motion of a particle with total 
mechanical energy FE satisfies the inequality 


E — U(z) > 0. 
The end-points of this region, for which 
E — U(z) = 0, 


are called the turning points of the motion. 


Section 4: The fabulous perfect spring 


For a system such as a pendulum or a toy car on a curved 

track, the law of conservation of energy has the form 
zmq + U(q) = £, 

where q is the displacement of the pendulum bob along 

its path or of the car measured along the track. For 

small oscillations about the lowest point of the path, 

q = qo, the potential energy function can be 

approximated by the first three terms of its Taylor series. 

This approximation leads to the prediction of simple 


harmonic motion of angular frequency w ~ ,/U’(qo)/m 


and period Tt ~ 27,/m/U"'(qo). 


Unit 8 Damped and forced vibrations 
Section 1: Damped vibrations 


1. Damping car be modelled by a damping force of 
magnitude r|z|, where r is a positive constant and z is 
the velocity of the particle, and direction opposed to the 
velocity of the particle. The SI units of the damping 
constant are Nm7's (or kgs‘). 


2. The equation of motion of a particle of mass m 
moving under the influence of a perfect spring of stiffness 
k and linear damping with damping constant r is 


mz+rz+kxr=0, 
where z is the displacement of the particle from its 


equilibrium position. Such a system is called a damped 
harmonic oscillator. 


3. For weak damping (r? < 4mk) the general solution 
of the equation of motion is 


a(t) = Ae” cos (Nt + ¢), 


where 
r V4mk — r? 
po ene SS. 
2m 2m 


All solutions of this form may be obtained with the 
arbitrary constants A and ¢ restricted so that A is 
non-negative and ¢ lies between —7 and 7. 

This motion is oscillatory with period tT = 27/Q. The 
amplitude of the oscillations decreases exponentially, so 
that 


ti + eity =e 


4. For critical damping (r? = 4mk) the general 
solution is 

a(t) = Be" +Cte™, 
where p = r/2m and B and C are arbitrary constants. 


This case gives the fastest return to the equilibrium 
position. 


5. For strong damping (r? > 4mk) the general solution 
is 
a(t) = Be "!* + Ce", 
where 
r— Vr? — 4mk 


2m 


r+~<Vr? —4mk 


and p, = a 


Pye = 


Section 2: Forced vibrations 


1. The equation of motion of a damped harmonic 

oscillator subject to an applied sinusoidal force of 

amplitude P and angular frequency w can be written as 
mizt+re+kx = Pcosut. 


The steady-state forced vibrations of this system can 
be found by using the phasor method. 


2. This equation of motion can be rewritten in the form 
: 2 os 
r+ 2awor + wor = — cosut, 
m 


where the undamped angular frequency w is defined 
(as in Unit 7) by 
k 
Wo = oye 
m 
and the damping ratio a is defined by 
r 


2/mk- 


= 


3. The steady-state vibration of the forced and damped 
harmonic oscillator is given by 


x = Acos(wt + ¢), 


where 
P/m 
A= ee si ee 
(we — w?)? + 4aPwew? 
and 
ae ae 
@ = — arccos 


(we — w?)? + 4da2wew? 


4. For a < 1/2, the amplitude of the forced vibrations 
exhibits a maximum at a particular value of the forcing 
angular frequency w. This effect is known as resonance. 


Section 3: The perfect dashpot 


1. A perfect dashpot represents a resistive force which 
is proportional to the relative velocity between two 
components of a mechanical system. 


2. The dashpot force has magnitude R = r|l|, where | is 
the length of the dashpot and r is a positive constant 
known as the dashpot constant. If the dashpot’s length 
is increasing, then the dashpot force is directed towards 
the centre of the dashpot; if the length is decreasing, then 
the force is directed away from the centre of the dashpot. 


3. The equation of motion of a particle attached to a 
perfect dashpot can be derived by assuming that the 
dashpot’s length is increasing. This equation of motion 
applies also to the case where the dashpot is contracting. 


Unit 9 Simultaneous linear algebraic 
equations 


In this unit the methods and theory are described for the 
special case of 3 equations in 3 unknowns, but the 
descriptions can easily be extended to the case of n 
equations in n unknowns. 


A set of 3 simultaneous linear algebraic equations 
in the 3 unknowns 21, £2, r3 is one that can be written 
in the form 
a2; + Opes ets = di, Ey 
a2@1 + bex2 + c2%3 = do, E2 
a321 + b3x2 + c3x3 = d3. E3 
The array of numbers 
al by CI: 
a2 bo €92 
a3 b3 G3 
is called the matrix (of coefficients) of the set of 
equations, and the array 
a1 by C1 dy 
a2 bo C2 do 
a3 bs C3 d3 


is called the augmented matrix. 


Section 1: Gaussian elimination 


A systematic procedure for solving such sets of equations 
is Gaussian elimination. 


Stage 1: Elimination 


(i) Eliminate x; from all but equation F1, by 
subtracting suitable multiples of FE, from each of 
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the equations below it. In the above example the 
multipliers used in calculating these multiples of Ey 
are ad2/a; and a3/a; for equations E2 and E3 
respectively. 


(ii) Eliminate x2 from all but the first two of the new set 
of equations by subtracting a suitable multiple of the 
new second equation from the equation below it. 


(iii) The new set of equations now has the ‘upper 
triangular’ form 


aay + bres + Crs Ay, Fy 
* * * 
b5 2X2 + Co%3 = ds, Eo 
* OK 
C373 = d3. E3, 


Stage 2: Back-substitution 
(i) Solve E3, for z3. 


(ii) Substitute this value of rz into E2,, and solve the 
resulting equation for x2. 


(iii) Similarly, substitute these values of x3 and z2 into 
FE, and hence find 7. 


Section 2: Types of solution 


1. A set of linear simultaneous equations will have: 
(i) a unique solution, or 
(ii) no solution, or 


(iii) an infinite number of solutions. 


To find out which case applies, look at the augmented 
matrix of the set of equations produced by stage 1 of 
Gaussian elimination. [In the table below, a matrix entry 
which must be non-zero is denoted by V, while a matrix 
entry which could be zero or non-zero is denoted by a 
star (*).] 


Characteristics of matrix Type of solution 
obtained after elimination 


Sf 
0 0 


No zero row of new 
coefficient matrix 


Unique solution 


(case (i)) 


Pao eet 
ican oe Se 
= See Bes 


A zero row of new coefficient 
matrix, but not of new 
augmented matrix 

(i.e. Equation F3, has form 
0 = non-zero number) 


No solution 


(equations inconsistent; 
case (ii)) 


* * 
See ie on ee 
Ga-2 4 


A zero row of the new 
augmented matrix 

(i.e. Equation £3, has form 
G= 0) 


An infinite number of 
solutions 


(case (iii)) 


In case (iii) we can obtain the solution set by taking 
x3 =k, where k is an arbitrary number, then 
back-substituting in Equations F2, and EF as usual. 
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2. The rows Ri, Ro,..., Rn of a matrix are said to be 
linearly dependent if some linear combination of 
these rows is zero, i.e. if constants ki, ke,...,kn (not all 
zero) exist such that 


k,R; + k2Re +--+ +krzRn = 0, 
where O here means a row of zeros. 


If the rows of the matrix A of coefficients in the given set 
of equations are linearly independent, then the solution is 
unique (case (i)). If the rows of A are linearly dependent 
but those of the augmented matrix A|b are linearly 
independent, then there is no solution (case (ii)). If the 
rows of A|b (and hence also those of A) are linearly 
dependent, then there is an infinite number of solutions 
(case (iii)). 


Section 3: Some possible difficulties 


1. (i) The kth pivot (or pivot element) is the one by 
which we divide to obtain the multipliers in the kth 
step of the elimination. 


(ii) The pivot row is the row in which the current 
pivot lies. 


(iii) Essential row interchange is the process of 
interchanging the rows of a matrix when the 
obvious pivot is zero. This can be done by 
interchanging the current pivot row with whichever 
of the rows beneath it will provide the pivot with 
the largest modulus. 


(iv) Partial pivoting is the process described in the 
preceding sentence with the difference that it is 
done even if the obvious pivot is not zero, provided 
that it will replace the obvious pivot by a pivot 
with a larger modulus. It ensures that no multiplier 
has modulus exceeding 1. 


2. The section shows some of the difficulties presented by 
simultaneous equations and their solution. 


(i) Ill-conditioning. A set of equations is said to be 
ill-conditioned if a ‘small’ change in the data (the 
coefficients, or the right-hand sides of the 
equations) produces ‘large’ changes in the solution. 
Ill-conditioning is a property of the equations 
themselves, not of the method used to solve them. 
One way to detect ill-conditioning is to vary the 
coefficients slightly and see if there is much change 
in the solution. A possible cure for ill-conditioning 

is to replace one of the equations by an equation 
containing some new independent information, and 
hence obtain a set of equations which are not 
ill-conditioned. 


(ii) Induced instability. This occurs when we start 
with a set of equations which are not ill-conditioned 
and yet the method used to solve the equations 
induces ill-conditioning in the problem during the 
process of solution. It is usually heralded by the 
appearance of a very large multiplier. Partial 
pivoting is used as a way of trying to avoid this 
(and to avoid the consequent build-up of error). 


Fe 


Section 4: The computer package 


Section 5: Special cases 


If the right-hand sides of the equations are all zero, the 
equations are said to be homogeneous. If the number of 
unknowns is the same as the number of equations, and 
the equations are linearly independent, then the only 
solution is the trivial solution, in which all unknowns 
equal zero. If, however, the equations are linearly 
dependent, then there is an infinite number of solutions. 


A set of m equations in n unknowns has an infinite 
number of solutions if m < n (unless the equations are 
inconsistent) and no solutions at all if m > n (unless at 
least m — n of the equations are linearly dependent on 
the others). 


Unit 12 Heat transfer 


Section 1: The nature of heat transfer 


1. Heat transfer means energy transfer by heating 
and/or cooling. This kind of energy transfer requires the 
existence of a temperature difference, and the energy is 
transferred from regions or bodies of higher temperature 
to those of lower temperature. There are three modes of 
heat transfer: conduction within a (normally stationary) 
material, radiation (or electromagnetic radiation) which 
does not require any material, and convection by 
motion of a fluid. In free convection the fluid motions, 
or convection currents, are driven by density changes 
caused by the temperature difference itself; in forced 
convection, they are driven by some other agent. 


2. Part of the energy of a body depends on its 
temperature. This part is known as its interna! energy, 
and the change of internal energy due to a change from a 
temperature 0; to a temperature 02 is mc(62 — 61), where 

m = mass of body, 

c = specific energy capacity 

or specific heat of the material. 

Power is the rate at which energy (of some particular 
kind) is produced. 


Section 2: Steady-state conduction in one 
dimension 


1. Fourier’s Law for steady-state conduction in one 
dimension is 


dé 
= —KA— 
q dz ? 
where 
q = heat transfer rate in the positive x direction, 


k = thermal conductivity of the material, 


A = area at right angles to the direction of heat transfer, 


9 = temperature, 
x = distance measured in direction of temperature 
variation. 


6 ° e 
— is known as the temperature gradient. 


dx 

Steady-state means that the temperature at any given 
point does not change with time; one-dimensional 
means that it depends on only one of the three possible 
space coordinates (in a suitable coordinate system). 


2. For steady-state conduction through a uniform slab of 
constant cross-sectional area A, Fourier’s Law reduces to 
0, — 02 

b 2 
where 6; — 62 (for 6; > 62) is the temperature drop 
across thickness b of the material. 


q=kA 


For steady-state conduction through a pipe of circular 
cross-section with internal radius r;, external radius r2 
and length /, Fourier’s Law leads to 
= 27K (44 as 62) 

1 “Tog. (r2/r1) 
where 6; — 62 (for 6; > 02) is the temperature drop 
across the pipe wall. 


Section 3: Walls and radiation 


1. Heat transfer through walls and roofs is often 
represented by steady-state one-dimensional models. 
The U value of a wall is defined as follows: 
rate of heat transfer per unit area 
U value = ——— 
overall temperature difference 
New techniques and components are being developed to 
satisfy the demand for very low U values in buildings. 


2. The thermal energy emitted by radiation from a body 
depends on the fourth power of the absolute 
temperature of the body, i.e. 

thermal energy emitted by radiation = 6%, 
where 

pu is a constant, 

64 = 273.2 + 6 = absolute temperature (in kelvins), 


6 = temperature in °C. 


Section 4: Convection and insulation 


1. Convective heat transfer in a fluid can be modelled 
using the following assumptions: 


(i) the temperature in the bulk fluid is uniform, say 6;; 


(ii) if the fluid meets a plane solid surface with area A 
and temperature @,, then the rate of heat transfer 
from surface to fluid is hA(@, — 07), where h is the 
convective heat transfer coefficient. 


2. The equation for heat transfer through a wall or 
window with area A is 

ee AU (6in 2 Pree 
where 6in and Oout (for Ain > Bout) are the bulk air 
temperatures inside and outside the building. Overall 
U values which take both conduction and convection into 


account can be calculated by adding formulas for 
adjacent temperature differences. 


In particular, for a solid wall 
eae 2 oe 
ae oe 
where 6 and « are the thickness and thermal conductivity 
of the wall and hin and Aout are the convective heat 
transfer coefficients at the inside and outside surfaces of 
the wall respectively. 
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3. Convective heat transfer from a surface such as the 
wall of a pipe can be increased by means of fins. The 
heat transfer in the fin can be modelled by assuming that 
the temperature in the fin depends only on its distance x 
from the pipe wall, dividing the fin into elements, and 
then considering convection from each element of the fin 
to the surrounding air and conduction from each element 
to its neighbours. This leads to the differential equation 
d*6 
dx? 
where h is the convective heat transfer coefficient, « is 
the thermal conductivity of the material of the fin, b is 
the fin’s thickness and 6, is the temperature of the air. 


2h 
= oe — 0.,), 


4. The heat transfer properties of a layer of fluid (such as 
the air space in double glazing) are complicated, and are 
best modelled using a combined heat transfer 
coefficient h., which is defined analogously to a U value 
but depends on the temperature difference across the 
layer as well as on its thickness. 


Unit 14 Vector algebra 
Section 1: Vectors: scaling and addition 


1. A (non-zero) vector is a mathematical object 
consisting of two parts: 


(i) a positive real number called the modulus or 
magnitude of the vector, and 


(ii) a direction in space (or in the plane). 


The modulus of a vector v is denoted by |v| or sometimes 
by uv. We also define the zero vector to have modulus 
zero; no direction is defined for it. The zero vector is 
denoted by O. In the context of vectors, we use the word 
scalar to denote a real number. 


2. Physical quantities that can be modelled by vectors or 
by scalars are called vector quantities or scalar 
quantities respectively. Examples of vector quantities 
include displacements (i.e. changes in position), velocities 
and forces. The displacement from a point P to a point 


Q is represented by the displacement vector PQ. 


3. Equal vectors have the same magnitude and direction. 
Scaling of a vector is multiplying it by a positive 
number; this leaves the direction unaltered and multiplies 
the magnitude by that number. Multiplying any vector 
by the number zero gives the zero vector. Multiplying a 
vector v by a negative number —m gives a vector —mv 
with magnitude m|v| and direction opposite to that of v. 
This new direction has the same orientation (i.e. in two 
dimensions, the same slope) as that of v but its sense is 
opposite to that of v. In particular the vector —v has the 
same magnitude as v and the opposite direction. 


4. A unit vector is one whose magnitude is 1. The 


, es: i ; ; 
vector ac is a unit vector with the same direction as a. 
a 


Unit vectors along the Cartesian axes are called 
Cartesian unit vectors, and are denoted by i, j, k. 
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5. Vector addition is defined by a parallelogram rule 


a 


or alternatively by a triangle rule 


a 


Vector subtraction is defined by 
a—b=a+(-b). 


6. Vector addition is commutative and associative (see 
Subsection 3.6 of this Handbook), and distributive with 
respect to multiplication by a scalar: 


(mi + m2)a = mia+t maa, 
m(a+ b) = ma-+ mb. 


Multiplication by a scalar has the ‘associative’ property 
that 


mi(mz2a) = (mimz2)a. 


The zero vector has the properties 0a = 0 and O0+a=a. 


Section 2: Components of a vector 


1. Given a Cartesian coordinate system, any 
two-dimensional vector a can be written uniquely in the 
component form 

a=aji+a2j, 
where i and j are unit vectors along the positive z- and 
y-axes respectively. The vectors aii and aaj are called the 
(Cartesian) component vectors of a, and the numbers 
a, and a2 are called the (Cartesian) components of a 
in the directions of i and j. The magnitude of a is 

a? + a2. To add two vectors, we add the respective 

components; to multiply a vector by a scalar, we multiply 
the respective components by that scalar. 


2. For three-dimensional vectors the representation 
a=aji+ a2j + a3k 


has analogous properties, where k is a unit vector along 
the positive z-axis, chosen at right angles to the x- and 
y-axes. A right-handed coordinate system is one 
satisfying the right-hand rule: if the thumb of your 
right hand points along i and the forefinger along Jj, then 
the middle finger can comfortably point along k. A 
right-handed coordinate system also satisfies the 
right-hand screw rule that a screw rotated from i to j 
will advance in the direction of k. 


3. If a point P has coordinates z, y, z relative to a given 
Cartesian coordinate system with origin O, the 


displacement vector OP is called the position vector of 
P, and is given by 
r= zi+yj+ zk. 


Section 3: Products of vectors 


1. The dot product or scalar product of two vectors 
a and b is the scalar defined by 


a+b = |al |b| cos 8, 
where 6 (for 0 < @ < 7) is the angle between the 
directions of a and b. 
In particular we have 

a-a=|al’. 
The angle 6 between two non-zero vectors a and b is 
given by 


2. The projection of a vector b onto a given direction is 
defined to be |b| cos 6, where @ is the angle between the 
direction of b and the given direction; it is equal to 
a-b/|a|, where a is any vector whose direction is the 
given direction (it is often convenient to make a a unit 
vector). Hence 


a-b = |a| times the projection of b onto the direction of a 


= |b| times the projection of a onto the direction of b. 


3. The dot product is commutative, 

a-b=b-a, 
and distributive, 

a-(b+c)=a-b+a-c, 

(at+b)-c=a-c+b°<c, 
but it is non-associative and there is no inverse. It 
combines with scaling according to the rule 

(ma) -b = m(a-b) =a: (mb). 
If a and b are non-zero vectors, then a- b = 0 if and only 
if a and b are perpendicular vectors. 


4. Dot products of Cartesian unit vectors are 
te eo ek = 1, 
Baa kek kt 0. 
If two vectors a and b have Cartesian components 
a1,a2,a3 and bj, b2,b3 respectively, then 
a- b= a,b; + a2b2 + azb3. 


The formula |a|? = a-a=aj+ az + a% is a particular 
case of this result. 


5. The cross product or vector product of two 
vectors a and b is the vector defined by 


a x b = |a| |b|sin@c, 
where 6 (for 0 < 6 < 7) is the angle between the 
directions of a and b, and c is a unit vector at right 
angles to both a and b. The sense of c is given by the 


right-hand screw rule, that a screw rotated from a to 
b advances along the direction of c. 


6. The cross product is distributive, 
ax(b+c)=axb+axece, 
(a+b)xc=axc+bxe, 

but is non-commutative, obeying instead the law 
a x b= —(b x a). 


It is non-associative, and there is no inverse. It combines 
with scaling according to the rule 

(ma) x b = m(a x b) = ax (mb). 
If a and b are non-zero vectors, then a x b = O if and 
only if a and b have the same orientation (i.e. they are 
parallel or anti-parallel). In particular, we have 
axa=0O. 


7. Cross products of the Cartesian unit vectors are 
ixteye tek xk = 0, 
Be ee 


Pe x 


If two vectors a and b have Cartesian components 
a1,42,a3 and bj, be, b3 respectively, then 


ab = (a2bs — a3b2)i + (az3bi — a1b3)j + (aibe2 — a2b; )k. 


8. The triple scalar product of three vectors a, b, c is 


a-(b xc). 


9. Given a plane, any vector n perpendicular to it is 
called a normal vector, and the equation of the plane 
can be written 


(2 — a)ni + (y — b6)n2 + (z —c)n3 = 0, 
where (a, b,c) is a point in the plane and 
n = n1i+ n2j + nz3k. If a and b are vectors lying in this 
plane, with different orientations, then a x b is a normal 
vector. 


Section 4: Differentiating vectors 


1. A vector function is a function whose domain is the 
set of all real numbers and whose codomain is the set of 
all vectors: to each real number t the vector function f 
assigns a vector f(t). Alternatively, the formula for f(t) 
may itself be described as a (vector) function of t. The 
Cartesian components of a vector function are real 
functions of the usual kind: 


f(t) = fi(t)it fo(t)j+ fs(t)k. 
The magnitude and direction of f(t) are also functions 
of t. A vector whose Cartesian components are constants 
(and hence whose magnitude and direction are constants) 
is called a constant vector. 


2. As t varies, the point whose position vector (relative 
to a given coordinate system) is f(t) traces out a curve in 
three dimensions, which helps in visualizing the vector 
function f. 


3. The derivative of a vector function f is 


(a (fie a ne | 


dt. hese h 
If f(t) has the Cartesian component form 
f(t) = fi(t)it+ fo(t)i+ fs(O)k, 
then 


df(t aitt).. dfett)—- tiatt 
in i+ : ‘i+ k. 
The rules for differentiating sums and products are 
d aw ._dv 
Fras +v) = eS 
d dm du 
a! u) = Bs + de’ 
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4 dt ps pigs 
= ak 


(ux v) =(S xv) + (ux =) 
dt ae qf 


4. The tangent at any given point on the curve 

representing a vector function f is in the direction of the 

df(t) 
dt 


a particle at time t, then its velocity vector is 


vector at that point. If r(t) is the position vector of 


dr(t) 
a If a 


vector function g(t) has constant magnitude, then 


dg 
Sie. See 
ea 


Unit 15 Newtonian mechanics in 
three dimensions 


Section 1: The ingredients of Newton’s second 
law 


1. In three-dimensional Newtonian mechanics the 
position of a particle is represented by the position 
vector 


r= 2i+ yj+ zk, 
where z, y and z are the coordinates of the particle, and 
i, j and k are unit vectors along the axes of a static 


right-handed Cartesian coordinate system that is 
calibrated in metres. 


2. The velocity vector v is found by differentiating the 
position vector with respect to t, the time measured in 
seconds: 

v=r=z2i+ yj+ zk. 
The magnitude of the velocity vector is the speed of the 
particle. The direction of the velocity vector is the 
direction of motion of the particle. 


3. The acceleration vector a is found by 
differentiating the velocity vector with respect to t: 


a=V=f = Zi+ yj+ Zk. 
This describes the rate of change of velocity both in 
magnitude and direction. The acceleration will be 


non-zero if either the speed or the direction of motion 
changes. 


4. Each individual force acting on a particle is 
represented by a force vector. If the direction of the 
force makes angles a, § and y with the positive z-, y- and 
z-axes, and if the magnitude of the force is |F|, then the 
force vector is 


F = |F|(cosai+ cos (@j+cosyk). 


The cosines that appear in this equation can sometimes 
be simplified by using the trigonometric identities 


cos (5 - $) = sing, cos (F +4) = — sing, 


cos(7 — ¢) = —cos¢. 


5. The gravitational force vector for a particle of 


mass m is mge, where e is a unit vector pointing 


vertically downwards and g ~ 9.81ms °. 
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6. If a particle experiences simultaneously a number of 
individual forces, represented by the vectors Fi, Fo,..., it 
behaves just as if it were acted upon by the single force 


F=F,+Fe+:::. 


This formula is known as the law of addition of forces, 
and the vector sum F is known as the total force 
vector. The components of the total force vector are 
denoted by F,, Fy and F;,. 


Section 2: Newton’s second law 


1. Newton’s second law, in its vector formulation, 
states that, at each instant of time, the motion of a 


particle is governed by the equation 
*r 
mr=F, i.e. m—~ =F 


dt? 


or, in terms of components, 
me = F,,° wa fy, - m2 = 25, 
where 
m is the mass of the particle, 
r = £i1+ yj + Zk is the acceleration vector, 
F = F,i+ F,j+ F.k is the total force vector. 


2. A particle that moves in a circle with constant speed 
is said to perform uniform circular motion, and its 
acceleration is directed towards the centre of the circle. 
For example, a particle of mass m for which 

r(t) = I(coswti-+ sin wt j), 
where w is a constant, performs uniform circular motion 
at speed lw round a circle of radius |. Its acceleration is 
—w*r(t) and the force on it is -mw7r(t), which has 
magnitude mw”*l and is directed from the particle towards 
the centre of the circle. 


Section 3: Algebraic equations and sloping 
tables 


1. A static particle, i.e. one which is at rest over some 
interval of time, experiences zero total force. 


2. A particle in contact with a solid surface experiences a 
normal reaction force, Fy, and a force due to friction, Fy. 


3. The normal reaction rule: Fy is perpendicular to 
the surface and |F y| is just sufficient to prevent the 
particle from sinking into the surface. For the particle to 
remain in contact with the surface, |Fy| > 0. 


4. The force of friction, or frictional reaction, F;, is 
parallel to the surface, in the direction that opposes 
motion. For a static particle, |F;s| is just sufficient to 
prevent sliding over the surface, but the particle can 
remain static only so long as |F | < |Fn|, where p is the 
coefficient of static friction. For a moving particle, 
[F;| = p'|Fn|, where py’ is the coefficient of kinetic 
friction. 


Section 4: Uncoupled differential equations 
and shot-putters 


1. We can model a projectile, such as the shot, as a 
particle that experiences only the downward force of 
gravity. The path or trajectory of a moving point is the 


set of points it passes through. The path of the shot is a 
parabola, that is, the graph of a quadratic function 
(with the y-axis vertical). 


2. If the shot-putter launches from zero height, the 
maximum range is 
Umax 
9 
where Umax is the greatest speed with which he can 
release the shot. The maximum range is achieved for an 
angle of release of a1 radians. 
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3. If the shot-putter launches from a height h, the 
maximum range is 


LJ 1 +3h/E; 


where L = vixax/g, and this is achieved for an angle of 
release of 


1 


J1 4+ 2h/L 


arctan 


Section 5: Coupled differential equations and 
pendulums 


1. This section discusses a particle of mass m suspended 
from a fixed point by a light taut inextensible string, 
where ‘light’ means ‘of negligible mass’, ‘inextensible’ 
means that 1, the length of the string, is constant, and 
‘taut’ means that the distance between the two ends of 
the string is l. 


2. If the particle moves in a horizontal circle, then its 
speed is constant, the tension in the string is constant, 
and the equations of motion for x, y and z are uncoupled 
(i.e. the equation for < does not involve y or z, and so on). 


3. If the particle moves in a vertical circle, then its speed 
is not constant, the tension in the string is not constant 
and the equations of motion are coupled. In this case the 
angle 6 between the string and the downward vertical 
satisfies the equation of motion 


9 = a sin 0. 
This is consistent with the energy equation 
im(10)° + mgl(1 — cos@) = EF = constant. 
The tension in the string is 


f = mile + mg cos @ 
2E 
= mg(3cos0@ — 2) + cc 
so long as the right-hand side of this equation is positive. 
As soon as J’ = 0 the string goes slack, and from then 
onwards a different description of the motion must be 


used. 


4. The sum of the kinetic and gravitational potential 
energies of a particle is conserved provided the particle 
experiences only the force of gravity and forces (like the 
tension in a string, the normal reaction force or the 
magnetic force on a charged particle) that act at right 
angles to its velocity. 


Unit 17 The dynamics of 
many-particle systems 


Section 1: Two-particle string—pulley systems 
Section 2: Systems of particles 


1. Newton’s third law states that the inter-particle 
forces acting mutually between a pair of particles are 
equal in magnitude and opposite in direction. If the 
particles are labelled 7 and 7, then this law may be 
expressed in terms of vectors by the equation 

Pig = —F iy, 
where F;,; is the force on Particle 2 due to the effect of 
Particle 7, and vice versa. 


2. For a system of n particles with masses 
m1,™M2,...,Mn and respective position vectors 
r1,Y%2,...,I%n, the centre of mass has position vector 


R= (mini + mere +--:+ Marr) = ig 
where 
M=mit+m2+-:-+mn= >> mi 
4=1 


is the total mass of the system. If the particles are acted 
upon by external forces F§*', FS**,..., F&** respectively, 
and if the inter-particle forces satisfy Newton’s third law, 
then the equation of motion of the centre of mass is 


MR me Fe*t 


where 
Fext = ae es po Beit obs po = = pe 
t=1 
is the total external force acting on the system. 


3. In the absence of external forces, the equation of 
motion of the centre of mass reduces to R = 0, indicating 
a uniform motion 


R(t) = Vot+ Ro, 
where the constant vectors Vo and Ro are respectively 


the velocity of the centre of mass and its initial position 
(at time t = 0). 


4. The centre of mass of a symmetric homogeneous 
object is located at the geometric centre of the object. 
The centre of mass of a collection of objects can be found 
by treating each object as if it were a single particle 
situated at its own centre of mass. 


Section 3: Linear momentum and energy 


1. The linear momentum p of a particle with mass m 
and position vector r is defined by 


D = ir: 
In terms of linear momentum, Newton’s second law 
can be expressed as 

p=F, 


where F is the total force acting upon the particle. In the 
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absence of any force on the particle, its linear momentum 
is conserved (that is, p is constant). If the mass m of the 
particle is constant, then Newton’s second law reduces to 


mr =. 


2. For a many-particle system, the total linear 
momentum P is the sum of the momenta of the 
individual particles. It is related to the total mass M of 
the system and the centre of mass R by the equation 


P = MR. 
Newton’s second law leads to the result 
P Ss pext 
= ’ 
where F°** is the total external force acting on the 
system. In the absence of external forces, the total linear 
momentum is conserved (that is, P is constant). If the 


total mass M of the system is constant, then the last 
equation reduces to 


MR = F**. 


3. For a system of particles involved in a collision in the 
absence of external forces, momentum conservation 
implies that Ppefore = Pafter. In the two-particle case, 
this is equivalent to 


mi1Ui + M2U2 = M1V1 + Me2V2, 


where the particles have masses m1, mz, initial velocities 
Uj, U2 and final velocities vi, v2, respectively. 


4. If a two-particle collision is elastic, then the kinetic 
energy equation 


3m |ui|’ + 3m2|u2|* = 5milvil” + 3malvel” 


is also satisfied. Otherwise the collision is inelastic. 


5. Newton’s law of restitution for a two-particle 
collision is 

(v2 —vi)-1= —e(u2 — u;)-], 
where | is a unit vector in the direction of individual 
particle momentum change, and the constant e is the 
coefficient of restitution for the objects represented by 
the particles. For an elastic collision, the coefficient of 
restitution is e = 1. If the collision is one-dimensional, 
then this law reduces to 


v2 — V1 = —e(u2 — U1). 


Section 4: Newton’s third law 
Section 5: Rocket motion 


1. The velocity v of a rocket satisfies the rocket 
equation 

ny =n 
where the rocket (with unspent fuel) has mass m, is acted 


upon by an external force F®*' and emits exhaust gases 
with relative velocity u. 


2. In the absence of any external force, the rocket 
equation reduces to 
m 


v= —u. 
m 
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Unit 18 Polynomial approximations 


Section 1: Basic theory 


1. A polynomial of degree n is a function that can be 


expressed in the form 
n n—1 
O,2 -+-On=1f Bei oe! 5 ee ame? 
where a, #0. The constant numbers @n,dn—-1,..., 41,40 
are called the coefficients. 


2. To evaluate the above polynomial for a given value of 
x by nested multiplication, use the recurrence relation 


yy = Tye Ene ee 8 12, ee 


with uo = an. Then up is the value of the polynomial. 


3. For a given function f(x) and a given number a in its 
domain, the nth Taylor polynomial for f about a is 


f(a) + (2-0) f'(a) + 5 (ea)? f"(@) 
1 n 
eae =e 


4. Taylor’s Theorem. Let the function f and its first 
n +1 derivatives exist and be continuous at all points 
between and including two given numbers a and z. Then 
we have 


f(x) = p(x) + (2), 


where 


a)* f (a). 


ple) = > ae -a)* f(a) 
k=0 
is the nth Taylor polynomial for f about a, and r(x), the 
remainder, satisfies 
Ce) = eae pe a ee 
for some number c, between a and z. Taylor’s Theorem 
with n = 0 is also called the Mean Value Theorem. 


5. The error function for the approximation of a given 
function f(x) by a polynomial p(z) is 


e(x) = p(x) — f(a). 
In the case of Taylor approximations it is related to the 
remainder by 


e(x) = —r(z). 


6. Suppose f is a given function and a,b, M are numbers 
such that 


Peps 
for all c satisfying a <c <b. Then, for any a satisfying 
a <a < ), the error in the nth Taylor polynomial 
approximation to f about a satisfies the error bound 
formula 


eh ie= al"t*'M fora<2<b. 


1 
(n+1)! 
Section 2: Two applications of Taylor 


polynomials 


1. A root of the equation f(x) = 0 means a solution of 
that equation. 


2. To calculate a root of the equation f(x) = 0 by the 
Newton—Raphson method, use the recurrence relation 
f(r) 

Lrtii =r f' (ar) . 
If the initial guess xo is not too far out, the sequence 
L0,21,22,... will converge to the root = e. for sufficiently 
large r, the value of z, is arbitrarily close to the root). 
The error in x,+1 is proportional to the square of the 
error in X,; so if the sequence converges at all it converges 
rapidly. 


3. The Taylor series method of order n for the 
solution of the first-order differential equation 

y' = m(z,y), 
where the value of y(xo) is given, involves the following 


steps. 


(i) Differentiate the expression m(z,y) with respect to 
x (remembering that y is a function of z) n — 1 
times to obtain expressions for 


dm d?m d"-!m 
Gn tr ¥)s dp? Pore i, my Guna (24): 


(ii) Choose a step size h. 


(iii) Calculate Yi, Y2,... from the recurrence relation 


Ypa1 = Y¥> + hY, + shy,” forest ShY,”), 
where : 3 

Lr = fo + Th, 

ee oe oe 


gti 23. 
= De, Ys) 


= dr”-1 


with the initial condition 
Yo = y(xo). 


(iv) Y;, is an approximation to y(z;,). 


Section 3: The catenary 


1. Hyperbolic cosines and sines are functions defined 
by 

coshae-5fe> + ~~); 

sinhz = $(e” —e ”). 
They have the properties 


2 ee 
cosh* x — sinh* z = 1, 


d 
—coshz = sinha, 


dx 

a= 

— sinhxz = cosha, 
dx 


2 4 6 


‘ 2s 4b 6 fb 

ee or a gee 
. x x? gt 

sin Cae ta a 


These functions are related to the hyperbola x? — y” = 1, 
in the same way that the ordinary cosine and sine 
functions are related to the circle. 
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2. The length L of a curve with equation y = y(z), 6. Simpson’s method for integrals, based on 
between two points on it with x-coordinates a and J, is approximating the area under a segment of a curve by a 
b 2 parabola, is 
a / 1+ ($*) dz. b 
. hae [she ae = $np(e0) + 44 er) + 2fle0) + 49( 

3. The curve with equation +2f(v4)+---+4f(an_-1) + f(rn)), 

= z cosh az, where 

: h = (b—a)/n, 


where a is any positive constant, is called a catenary. It 


models the shape of a rope suspended between two points. tp=a+rh (r=0,1,...,n), 


and n must be even. 


Section 4: Interpolation and integration 7. The three integration formulas are examples of 
Newton—Cotes formulas, obtained by integrating an 

1. If the values of a function f(x) are known at points interpolating polynomial. They are said to be in 

L0,11,..-,Xn, arranged in increasing order but not composite form if n > 1 for Euler’s and the trapezoidal 

necessarily equally spaced, then the use of some other formulas, and if n > 2 for Simpson’s. For a given step 

function g(x), which takes the same values as f(x) for size h, Euler’s formula is normally the least accurate of 

L = 2X0,L1,..-,Ln, to approximate f(z) for values of x the three and Simpson’s the most accurate. 


between 2p and Zn, is called interpolation; the 
corresponding procedure for values of x less than xo or 


greater than x, 1s called extrapolation. 
Unit 19 Numerical methods for 


2. There is a unique interpolating polynomial of differential equations 
degree < n which takes the same values as f(z) for 
T= £0,%1,.--, En. ltsn+1 coefficients, =. In all these methods the differential equation is taken to 
Qo, @i1,..-,@n, can be found by solving the n + 1 — 
simultaneous linear equations ; 
ao +aixi+---+anz; =f(ai) (i=0,1 n) Pe oe 
ee = : eee with y(zo) given, and we seek approximations to y(z) for 
¢-= 21, 05,.. . where 


3. Provided xo < x < an, the error function ¢(x) for the 

approximation of f(x) by the interpolating polynomial 

defined above satisfies and h is the step length or step size. The approximation 
1 


fy = Xo + rh 


As\ == fe a) (a — 21) rege he which means y(z,); is denoted by Yr. The 
(n+ 1)! approximation to y,, which means y’(z,), is denoted by 
where cz is some number satisfying ro < cz < Xn. If M is Y, and is given by the formula 
a number such that Y,/ = m(zr, Y>). 


oa 


for all c between zo and zp, then ¢(z) satisfies the error 
bound formula 


1 
le(z)| < Ga 


Section 1: Numerical methods 


1. By an integration method for the numerical 
solution of differential equations we mean a method 
based on a numerical integration formula. 


(xv — zo)(x@ — 41)... (x@ — 2n)|M. 


4. Euler’s integration formula, based on 


approximating the area under a segment of a curve by a 2. Euler’s method uses the recurrence relation 
rectangle, is Gus th. 
b 
: Jiaide = hf a0) ee ee 3. The trapezoidal method uses the formula 

where fee Y,.'+ +h(Y, + Y,44). 

h=(b—a)/n In the case of a linear differential equation, for which 
and 20,21,...,2n—1 are equally-spaced points defined by m(z,y) =lU(x)y+ k(2), 

o = 464 th = 0A. a the formula of the trapezoidal method can be solved for 

Y;-41 to give 
5. The trapezoidal method for integrals, based on . (1+ 5Al,)¥, + 5h(ky + krp+1) 
r+1 = y 


approximating the area under a segment of a curve by a 
trapezium, is 


b 
Z f(x) da ~ h(3 f(zo) + f(ai) +--+ + f(tn-1) + 3 f(2n)), Sf 
= 4. An explicit method is one where the recurrence 


where relation for Y;41 does not involve Y,,,, Y;41, or higher 

h = (b—a)/n, derivatives at z-41. An implicit nono is one where 
the recurrence relation does involve Y,,, or higher 
derivatives at %r41. 


where |, = I(z,), ky = k(az,), and so on. 


f= rh fr = 
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5. The local truncation error is the error introduced 
in a single step, estimated as the value of Yr41 — yr+1 
when Y; = y,. The principal term in the local 
truncation error is the first non-vanishing term of the 
Taylor series about xz, in powers of h for the local 
truncation error; this series is calculated by subtracting 


the Taylor series 
ae 2 


Yeti = Yr t hye + Gh ye + 
from the corresponding series for Y,+41 obtained from the 
recurrence relation. 


6. The principal term in the local truncation error for 
Euler’s method is —ih*y’, for the Taylor series method 


of order n (with n > 2) it is ant ty Hn +1) aed 
‘tt 


for the trapezoidal method it is ah? y: 


Section 2: The predictor-corrector method 


1. The Euler-trapezoidal predictor-corrector 
method calculates Y,+1 in two stages: 


re Heat a where Y, = m(z,r,Y,), 


Yri1 =Y¥,r + $R(¥;+Y741) where YY) = m(zr41, Yr41). 


The principal term in the local truncation error for this 
method involves h?. 


2. In general a predictor-corrector method is one 
where we calculate Y;11, the approximation to yr+1, in 
two stages: the first stage uses an explicit method to give 
a first approximation Y,",,, which is then used on the 
right-hand side in the recurrence relation of a more 
accurate implicit method to give an improved 
approximation Y;4+1. 


Section 3: The analysis of numerical methods 


1. A one-step numerical method for the solution of 
the differential equation y’ = m(z,y) is a method which 
approximates the differential equation by a first-order 
recurrence relation of the form 


Yr41 = Yr + h o(zr, p Sole fee h). 
It is consistent with the differential equation if 


P(Er, Yr, Yr, 0) = M(Lr, Yr). 
An alternative statement of the definition of consistency 
is 

‘ p Se StS Yy. 

lim ————— = pc ter 
h20 h es ) 
2. For a fixed value of z, say xz“, the.global error of a 
numerical method for calculating y(x*) using N steps of 
length h = (x* — xo)/N is defined as Yn — y(z"). 
The method is said to be convergent on some interval 
[xo, b] if for all z* in this interval (i.e. all 2” satisfying 
ro < x* <b) we have 


lim Yy = y(2”*). 
lim Yn = y(z") 


3. A one-step method is convergent if and only if it is 
consistent. Further, if the principal term in the local 
truncation error is h?*', for some integer p, then the 
global error at x” is given approximately by 


Yn — y(x") ox Ch’, 


where, for small h, C does not depend on h. 


Fete 


In particular, the global error for Euler’s method is 
approximately proportional to h, whereas the global 
errors for the trapezoidal method, the Euler-trapezoidal 
method and the Taylor series method of order 2 are all 
approximately proportional to h? for small h. 


4. A one-step method, applied to a given differential 
equation problem, is said to be absolutely unstable if 
the resulting recurrence relation problem is absolutely 
ill-conditioned, and to be absolutely stable if the 
recurrence relation problem is absolutely 
well-conditioned. 


5. We can test the stability of a given numerical method 
by applying it to the test problem y’ = ay with 

y(xo) = 1, where a is a constant whose value is to be 
chosen. The interval of absolute stability of the 
method is the set of values of ha for which the method is 
absolutely stable. 


6. The interval of absolute stability for Euler’s method, 
for the Euler-trapezoidal predictor-corrector method and 
for the Taylor series method of order 2 is (—2,0), i.e. 
these methods are stable if —2 < ha < 0. The interval of 
absolute stability for the trapezoidal method is (—oo, 0), 
i.e. the method is stable for all negative values of ha. 


7. Applied to the equation y’ = m(z, y), a given 
numerical method is absolutely stable if h——- (ae, Y,-) lies 


in the interval of absolute stability for all values of r used 
in the calculation. 


Section 4: Simpson’s method 


1. Simpson’s method uses the recurrence relation 
Yr41 = Yp-1 + $A(Y,_1 + 4Y; + Y,41)- 


Since the recurrence relation is of second order, this may 
be called a two-step method. 


2. Simpson’s method is implicit, but if it is applied to a 
linear differential equation 
y =I(x)y + k(x), 
the recurrence relation can be solved for Y;+1 to give 
Ahl,-Y> + (3 + Alr—1)¥p—1 + A(ke—a + 4ler + br+1) 
3 — Aly4i 


where I, = I(z,), k- = k(z,-), and so on. 


bre = 


3. The principal term in the local truncation error for 
Simpson’s method contains the factor h°, and the method 
is therefore capable of very high accuracy. 


4. A spurious solution of the recurrence relation for 
Simpson’s method is a term in the general solution which 
does not correspond to any solution of the differential 
equation. For the differential equation 


y =U(x)y + k(z), 
the spurious solution is harmless provided I(x) is positive, 


but if I(x) is negative the spurious solution makes the 
method unstable. 


Unit 20 Matrix algebra and 
determinants 


Section 1: Some simple matrix operations 


1. An m X n matrix is a rectangular array of numbers 
called elements, in m rows and n columns. The two 
numbers m and n together specify the size of the matrix. 
The element in the ith row and jth column of the matrix 
A can be denoted by a;;. A column matrix or column 
vector is a matrix with just one column; a row matrix 
or row vector has just one row. 


2. Two matrices A and B are equal if they are the same 
size and corresponding elements in each matrix are equal. 


3. If A is a matrix and k is a real number, then kA is 
the matrix of the same size whose elements are k times 
the corresponding elements of A. 


4. If two matrices A and B are the same size, then we 
can add them, the sum A+ B being the matrix formed 
by adding the corresponding elements of A and B. 


5. A matrix whose elements are all zero is called a zero 
matrix and is written 0. For any matrix A we have 


A+0=0+A=A 
and CA = 0. 


6. Given any matrix A, the matrix (—1)A, whose 
elements are the negatives of those of A, is written —A. 
This matrix has the property 


A +(—A) =0, 
and it can be used to define matrix subtraction as 


A-—-B=A+(-B). 


7. The algebraic operation of adding matrices of a given 
size is closed, commutative and associative, and is 
distributive with respect to multiplication by a scalar (see 
Section 3.5 of this Handbook). 


Multiplication by a scalar obeys the further rules 
(kiko)A = ki(k2A), 
(ky + k2)A —kA+koA. 


Section 2: Matrix multiplication 


1. If A is an m X n matrix and B is an n X p matrix, 
then we can multiply them, with A written on the left. 
The product C = AB is an m x p matrix whose element 
in the ith row and jth column is 


n 
Cj = ) dikbp;- 
p=t 


2. Matrix multiplication is associative, and is distributive 
with respect to matrix addition, i.e. 


A(BC) = (AB)C, 
A(B+C)=AB+AC 

and (A+B)C=AC+BC. 

However, it is not (in general) commutative: 


BA F AB. 
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3. The product of a matrix and a column vector, with 
the matrix on the left, is another column vector. A set of 
simultaneous linear algebraic equations can be written in 
the form 


Ax = BD, 


where A is the (given) matrix of coefficients, x is a 
column vector whose elements are the unknowns, and b is 
a given column vector. Moreover, if A is a given matrix 
and x is a variable column vector of suitable size, then we 
can say that A transforms x to Ax (i.e. A defines a 
function whose domain is a set of column vectors; this 
function is called a linear transformation). 


4. The transpose of a matrix A, written A’, is 
obtained by interchanging the rows and columns of A, 
i.e. the element in the ith row and jth column of A? is 
a;;. The main properties of transposed matrices are 


(ary = A, 
(Aj As A Se ACA: 


Section 3: Change of axes 


Suppose Oz, Oy and O’z’,O’y’ are two Cartesian 
coordinate systems in the plane, with origins O and O’, 
that the origin O’ of the second system has coordinates 
(a,b) with respect to the first system and that the 
direction of the O’z’ axis is obtained from that of Ox by 
an anticlockwise rotation through an angle a. Then the 
coordinates (x,y) of a point P with respect to the first 
system and its coordinates (x’, y’) with respect to the 
second are related by 


| |} cosa Sine | 42 —a 
y |  |—sina cosa||y—6]- 


Section 4: Square matrices and their inverses 


1. A square matrix is one with the same number of 
rows as columns. The main diagonal of a square matrix 
A consists of the elements ai1,@22,.... 


2. A unit (or identity) matrix, written I, is a square 
matrix with elements 1 on the main diagonal and 0 
everywhere else. For any square matrix A and unit 
matrix I of the same size, we have 


AI=IA=A. 


3. A singular square matrix is one whose rows are 
linearly dependent. A non-singular square matrix is one 
whose rows are linearly independent. Every non-singular 
square matrix A has a unique inverse A~’ with the 


property 

AAS AS Ave. 
Singular matrices (and non-square matrices) do not have 
inverses. 


4. To compute A’, start with the matrix pair A|I and 
perform the same Gaussian elimination type row 
operations on both matrices, including multiplying or 
dividing a row by a non-zero constant. The linear row 
operations which convert the left-hand matrix A to a unit 
matrix I will convert the right-hand matrix I to A~?. 
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In particular, the 2 x 2 matrix 
Bee 
has inverse 


Se 4S 
ie -aE|-2 di 


5. Matrix inverses have the property 
(Aj Ao Ap aA =A A A A. 


6. A square matrix A is said to be symmetric if 
A? =k 


Section 5: Introduction to determinants 


1. To every square matrix A there corresponds a number 
called its determinant, written det A or |A|. 


2. If A is a 2 X 2 matrix, then 


Q11 ai2 
Q21 a22 


det A = 


= 4211422 — 412421. 


3. If A is a 3 x 3 matrix, then 


Q11 @19° “@zs 
det A =]a21 a22 423 


= 20°: G95 21 G23 
ae eS == Ero 
Qa32 433 Qa31 a33 
Q21 a22 
+ @13 
Qa31 432 


An n X n determinant can be defined in terms of 
(n — 1) x (n — 1) determinants in a similar way. 


4, The rows (or columns) of a square matrix A are 
linearly dependent if and only if det A = 0. Therefore 
A~? exists if and only if det A # 0. 


5. Linear combinations of rows (or columns) of a 
determinant can be added to any other row (or column) 
without changing the value of the determinant. 


6. In an upper triangular matrix, all elements below 
the main diagonal are zero; in a lower triangular 
matrix, all elements above the main diagonal are zero; in 
a diagonal matrix all elements not on the main 
diagonal are zero. The determinant of an upper 
triangular, lower triangular or diagonal matrix is the 
product of the elements on its main diagonal. 


7. To evaluate a determinant, use the operations of 
Gaussian elimination, namely adding or subtracting a 
multiple of one row from another and including row 
interchanges if necessary, to reduce the matrix to upper 
triangular form. Provided no row interchanges have been 
made, the upper triangular matrix has the same 
determinant as the original matrix; every interchange of a 
pair of rows reverses the sign of the determinant. 


Unit 21 Eigenvalues and eigenvectors 


Section 1: The theoretical eigenvalue problem 


1. Given a square matrix A and a number \X, the 
equation 


Ax = A’x 


for the column vector x may have solutions other than 
the zero column vector. If it does, then we say that the 
number A is an eigenvalue of the matrix A, and that 
the non-zero column vector x is an eigenvector of A. 
Any non-zero scalar multiple of an eigenvector is also an 
eigenvector, so there are an infinite number of 
eigenvectors of A corresponding to each eigenvalue. 


2. One way of finding the eigenvalues A of a given square 
matrix A is to solve the characteristic equation 
det(A — AI) = 0, 
which is a polynomial equation in 4. Eigenvectors can 
then be found by solving the equation 
(A — AI)x = 0 


for x. 


3. If the characteristic equation has n different solutions 
A1,..+,;An, where the size of the matrix A is n x n, then 
the eigenvalues of A are said to be distinct. If not, then 
the polynomial det(A — AI) has at least one repeated 
linear factor (A — A;)” with m > 2, and the 
corresponding eigenvalue 4; is said to be repeated. 


4. If the eigenvalues of an n x n matrix are distinct, then 
the matrix has n linearly independent eigenvectors. If the 
eigenvalues are not distinct, this may or may not be the 
case. 


Section 2: Iterative methods for finding 
selected eigenvalues 


1. If the square matrix A has eigenvalues A1, A2,...,An 
then for any real p, q: 


(i) A-+gqI has eigenvalues 41 + q,...,An + q; 

(ii) if A is non-singular, then A~* has eigenvalues 
yee. ie 

(iii) if A — pl is non-singular, then (A — pI)~* has 
eigenvalues (Ai; — p)~*,...,(An —p) 3 

(iv) A, A+qI, A~‘ and (A — pI)’, if they exist, all 
have the same eigenvectors. 


2. Direct iteration. We use the recurrence relation 


Vraui = Ay,/Gr41, 
where a,+1 is the element of largest modulus in Ay,. If 
Y, is not an eigenvector and not the zero vector, then as 
r becomes large, a, approximates to the eigenvalue of A 
with largest modulus (assumed real and distinct) and y,, 
approximates to a corresponding eigenvector. 


3. Inverse iteration. We calculate y,,, either from 
its definition y,,, = A~‘y,, or by solving the 
simultaneous equations 


AY, 44 =r 


and then obtain y,,, from 


tego V1] Orsi, 
where a,+1 is the element of largest modulus in y,,,. 
For an arbitrary non-zero starting vector y, which is not 
an eigenvector, and with A non-singular, the number 
1/a, approximates to the eigenvalue with smallest 
modulus (assumed real and distinct) and y,. approximates 
to a corresponding eigenvector, for sufficiently large r. 


4. Modified inverse iteration. ‘To find the eigenvalue 
closest to p and a corresponding eigenvector of A, we 
start with an arbitrary non-zero vector y, which is not an 
eigenvector and calculate y,,, either from its definition 
Y-41, = (A -—pl)“*y,. or by solving 


(A — pl)y,41 = Ve 
and then obtain y,,, from 


bre ees Vp44/Ortis 
where a,+1 is the element of largest modulus in y/,,,. 
If A — pI is non-singular, and r is large enough, the 
number p+ (1/a,) approximates to the eigenvalue closest 
to p (assumed real and distinct), and y,. approximates to 
a corresponding eigenvector. 


5. For a 3 x 3 matrix where we know the eigenvalues of 
largest and smallest modulus, denoted here by Ai and 2, 
a good choice of p in modified inverse iteration is either 


5(\A1| + |A2|) or —5(|Ar| + |Aal)- 


Section 3: Decomposition methods for finding 
all the eigenvalues 


1. To decompose a square matrix A is to express it as 
a product of two other square matrices, and the LU 
decomposition of A (when it exists) is 

A=LU, 
where L is lower triangular with all diagonal elements 
equal to 1, and U is upper triangular. 


In particular, a 2 x 2 matrix A has an LU decomposition 


Qi. G12} _ 1 O} | air a2 
Q21 422 a21/a11 1 0 det A/ai ‘ 
provided that aii is not zero. 


x 


2. The eigenvalues of a diagonal, upper triangular or 
lower triangular matrix are equal to the elements on the 
main diagonal. 


3. If A, B and P are n x n matrices, with P 
non-singular, then: 


(i) AB and BA have the same eigenvalues; 


(ii) A and P~'AP have the same eigenvalues. Further, 
if A has n distinct eigenvalues \1,...,An, and P is 
the matrix whose columns are the corresponding 
eigenvectors of A, then 


xi 0 


P-'AP e 


() 


~ 
) 


4l 


where all off-diagonal elements of this last matrix 
are zero. Even if the eigenvalues are not all 
distinct, the result still holds provided a set of n 
linearly independent eigenvectors exists. 


4. The LR method for calculating eigenvalues. Given 


a square matrix A, form the sequence Ao, Ai,..., as 
follows: 

Ao =A, 

Ar+1 seas Le, 
where A, = L,U, is the LU decomposition of A,. For 
suitable A, the sequence Ao, Ai,... will converge to an 


upper triangular matrix whose diagonal elements are the 
eigenvalues of A arranged in order of decreasing modulus. 


5. If the LR method applied to a matrix A fails to 
converge, try applying it to the matrix A + qI, where q is 
a real number. The eigenvalues of A + gI are 

Ai +q,A2 + q,..., where Ai, A2,... are the eigenvalues 

of A. 


Unit 22 Simultaneous differential 
equations 


Section 1: First-order systems 


1. A linear first-order system of differential equations 
is one that can be written in the form 

Ai(t)x(¢) + Ao(t)x(t) = g(t), 
where x(t) is a column vector whose elements are the 
unknown functions of t, x(t) is the column vector 
obtained by differentiating these functions, Ai(t) and 
A2(t) are given matrices and g(t) is a given column 
vector whose entries may depend on t. 


A linear second-order system is one that can be 
written 


Ax (t)%(t) + Ao(t)x(t) + As(t)x(t) = g(t). 
If all the entries in the matrices A; and Ag, and As if 
present, are constants, we have a constant-coefficient 
system. If g(t) = 0 for all t, we have a homogeneous 
system. 


2. The normal form (if it exists) of a linear first-order 
constant-coefficient system is 


x(t) = Bx(t) + h(t), 


where B is a constant square matrix. 


3. If B is an n X n matrix with n linearly independent 
eigenvectors, say a1,...,@n, Corresponding to (not 
necessarily distinct) eigenvalues \i,...,An, then the 
general solution of the homogeneous system 


x(t) = Bx(t) 


x(t) = Cyaje™!’ +---+ Crane”, 


where Cj,...,Cn are arbitrary constants. If B is real 
then the expressions a,e*”* in the general solution are 
either real or occur in complex conjugate pairs. If ae 
and ae* is such a complex conjugate pair, then the real 
pair Re(ae*’) and Im(ae*’) may be used instead. 
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Section 2: Further methods for linear 
first-order systems 


1. If B is an n X n matrix with n linearly independent 
eigenvectors, then the system 


x = Bx + h(t) 
can be solved by putting 
x = Py, 


where P is a square matrix whose columns are the 
eigenvectors of B. The resulting equation for y is 


y =P 'BPy+P ‘h(t) 
and, since P~'BP is a diagonal matrix (its diagonal 
elements being the eigenvalues of B), we can solve the 


equations for the various elements of y individually. The 
solution x is then found from x = Py. 


2. Euler’s method. This numerical procedure can be 
used for first-order systems of the form 


x= m(x, tf), 
with initial condition x(to) = xo, where m is a column 


vector whose elements may depend on the elements of x 
and on t. 


The approximation to x(t,), where tr = to + rh with h 
the step length, is denoted by 


R= (Ais aie Read 
and we take Xo = Xo. The recurrence relation is 
Xt) =e So ae) 1 8) 
or X41 = X, + hme 7,): 


Section 3: Second-order homogeneous systems 


1. If B is an n X n real matrix with n linearly 
independent eigenvectors, the general solution of the 
system 


xX = Dx 
is a sum of the following terms: 


(i) for each real positive eigenvalue 4, with eigenvector 
a,, a term 
ar (Cre yet — De = 
(ii) for a zero eigenvalue with eigenvector a,, a term 
ar(C;, + Dt); 


(iii) for each real negative eigenvalue 4, with 
eigenvector a,, a term 


a,(C, cos \f —A-t + D; sin \/ —Art); 


(iv) for each complex pair of eigenvalues \, and \,, a 
term 


Ke Re(a,eV **) +C, Im(a;eV>") 
+D, Rete v=) + D, infec. 


where a, is an eigenvector for A-. 


Section 4: Forced oscillations 


1. Theorem 1: Suppose that x, is a particular solution 
of the system 


Aix + Aox + A3x = h(t). (1) 
Then the general solution of this system is 


X = Xp + Xe, 


where x, (the complementary function) is the general 
solution of the associated homogeneous system 


Aix + Aox + Asx = 0. 


2. Theorem 2: Suppose the matrices Ai, A2, A3 are 
real, constant, symmetric and square, and that all their 
eigenvalues are positive. (A symmetric matrix is one that 
equals its own transpose.) Then every solution of 


Aix(t) + A2x(t) + A3x(t) = J 
becomes small (i.e. every entry in the column vector x(t) 
becomes small) as t becomes large. 


3. To look for a particular solution of 
Aix + Aox + A3x = bi coswt + be sin wt, 
write the system in the phasor form 
Aix + Aox + A3x = Re(be*”’), 
where b = b, — ibg, and look for a solution of the form 
x = Re(ze*’), 
where z is a column vector whose entries are complex 
constants. 


Unit 24 Normal modes 


Section 1: Modelling vibrating systems 


1. A lumped-parameter model is a model of a real 
mechanical system in which each component of the model 
has only one of the three properties mass, stiffness and 
damping. The only components used in this unit are 
particles (having mass only) and perfect springs (having 
stiffness only). 


2. The number of degrees of freedom in a 
lumped-parameter model is the number of coordinates 
necessary to specify the configuration (i.e. the positions 
of all the particles). 


3. The behaviour of a vibrating system is said to be 
forced if time-dependent external forces act on it; 
otherwise the behaviour is said to be free. 


Section 2: Free undamped vibrations with two 
degrees of freedom 


1. The equilibrium configuration (or static 
configuration) of a lumped-parameter system is the 
configuration for which the total force on each particle is 
zero; an equilibrium position of a particle is its 
position when the system is in an equilibrium 
configuration. 


2. A normal mode of a mechanical system is a motion 
of the system in which all the particles execute simple 
harmonic motion with the same angular frequency. The 
angular frequencies for which this is possible are called 
normal mode angular frequencies. 


3. In normal mode motion, the displacements of any pair 
of particles from their equilibrium positions are in a 
constant ratio, called the displacement ratio for that 
pair of particles. If the displacement ratio is positive, the 
two particles are said to be in phase; if it is negative, 
they are said to be phase opposed. 


4. The equations of motion for an undamped vibrating 
system can, by dividing each equation by the relevant 
mass, be put in the form 

x = Hx, 
where x is a column vector whose entries are the 
displacements of the various particles from their 
equilibrium positions, and H is a square matrix with 
constant entries. Then each eigenvalue of H is —w”, 
where w is a normal mode angular frequency, and the 
corresponding eigenvector, if the top entry in it is chosen 
to be 1, gives the displacement ratios for the 
corresponding normal mode. 


5. The most general motion of a vibrating system is an 
arbitrary linear combination of the normal mode motions. 


Section 3: Vibration absorption 


1. Forced vibrations. When a sinusoidally varying 
external force is applied to some part of a 
lumped-parameter system, a motion is in general possible 
for which all particles execute simple harmonic motion at 
the angular frequency of the external force. This 
steady-state motion can be used to model the 
steady-state response of a real system to such a force. 


2. If the frequency of the external force is close to a 
normal mode frequency, then the forced vibrations 
exhibit resonance—the amplitude of the steady-state 
response is very large. 


3. There may be a frequency for which one of the 
particles remains at rest in the steady state. This 
phenomenon can be used to isolate the object modelled 
by that particle from the effects of forced vibration at (or 
near) this frequency. The system is then said to act asa 
vibration absorber. 


Section 4: Extending the scope 


A degenerate system is one for which one of the normal 
mode angular frequencies is zero. In the corresponding 
motion each particle moves with the same constant 
velocity (rather than sinusoidally). 


Unit 25 Functions of more than one 
variable 


Section 1: First-order partial derivatives 


1. Any expression or formula which involves two 
variables, say x and y, and whose value is uniquely 
determined by the values of x and y, is called a function 
of the two variables z and y. If a variable z is a 
function of x and y (i.e. z = f(x,y), where f is a function 
of the two variables x and y) then we call x and y the 
independent variables and z the dependent 
variable. The set of possible pairs of .values (x, y) for the 
two independent variables is called the domain of the 
function; it can be thought of as a region in the 


(x, y)-plane. 
Functions of three (or more) variables are defined 
analogously. 


2. A function of two variables, say f(x,y), can be 
represented in three-dimensional space by the set of 
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points whose coordinates z, y, z satisfy z = f(x,y). Such 
a set of points is an example of a surface. 


3. A section function is a function of one variable, 
obtained from a function of several variables by holding 
all but one of the independent variables fixed. 


4. If f(x,y) is a function of two variables, and a,b are 
two constants, then the (first-order) partial 
derivative of f(z,y) with respect to z at (a,b) is defined 
as the (ordinary) derivative at a of the section function 


f(z, b), and is denoted by (a, b). Similarly, the partial 


derivative of f(x,y) with respect to y at (a,b) is the 
derivative at b (with respect to y) of the section function 


f(a,y), and is denoted by Se (a, b). 


5. We often use the partial derivatives of f(x,y) at (z, y) 


instead of at (a,b). The partial derivative —(z, y) is 


calculated by differentiating f(z, y) with respect to z 
while treating y as a constant. This procedure is called 
differentiating f(z,y) partially with respect to z. 


Similarly, dy is calculated by differentiating f(z, y) 


partially with respect to y at constant z. 


First-order partial derivatives of functions of three or 
more variables are defined similarly. 


Section 2: Some applications of partial 
derivatives 


1. Let f(x,y) be a well-behaved function—that is, one 
such that the surface z = f(x,y) is smooth in the sense of 
having no folds or breaks. Then the Taylor polynomial 
of first order for f at (a,b) is 


pi(ey) = f(a.) + (2 — a) 52(a,8) + (y ~ 8) 52(a,) 


2. The first-order Taylor approximation formula at 
(a,b), f(x,y) ~ pi(z,y), holds for (x, y) close to (a,b). It 
can be written 

Of Of 


Flay) & F(a,b) + (@ ~ a) 5° (a,b) + (y 6) 5 (a,b) 


or 6F 22 oF bx + Sou 


where 


6f = Fiagy — fie), 
6x = xr —a, 


édy=y-—b. 


3. The triangle inequality is 
ja + b| < Ja] + [OI 


and holds for any two numbers a and b. 


4. The tangent plane to the surface z = f(x,y) at the 
point (a,b, f(a,b)) has the equation z = pi(z, y). 


5. The chain rule for the function of two variables 


f(x,y) is 
df _ Of de | af dy 
dt: Ox: dt: = dy-dt’ 


where zx and y are functions of the single variable t. 
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6. The above formulas can be generalized to functions of 
three or more variables. 


Section 3: Higher-order partial derivatives 
and Taylor polynomials 


1. The second-order partial derivatives of a function 
f(z, y) are 


Orr See 
se lay ae 
Of as Ofer 
a ee 


it ade 
OxOy Ox \ dy)’ 


Of Oe 
Oy2 Oy \ Oy) — 


Third- and higher-order partial derivatives are defined 
similarly, as are the partial derivatives of functions of 
more than two variables. 


2. The Mixed Derivative Theorem. If f(z,y) isa 
well-behaved function of two variables, then 


a’ f arf 


Ardy OyOxr 


3. The Taylor polynomial of second order for 
f(z,y) at (a,b) is 


pa(su) = f(a,b) + (2 ~ a) £(a,b) + (y =) F(a, 


2 0 , z. 
+ Hea) Fa lab) +(e - ally - Dy 
e? 
41g 6) 55 (a0) 


4. For (z,y) near to (a,b), the second-order Taylor 
approximation at (a,b) is f(z, y) ~ po(x, y) and is 
normally more accurate than the first-order Taylor 
approximation. 


Section 4: Maxima and minima 


1. If, for all x and y sufficiently close to a and B, 
f(x,y) < f(a, 5), 


then f has a (local) maximum at the point (a,b); but if 
the inequality is 


f(x,y) > f(a,6), 


then f has a (local) minimum at (a,b). 


2. An extremum is a point which is either a local 
maximum or a local minimum. 


3. A stationary point of a function f(z, y) is a point 
(a.b) where 


0 6) 
=F (a,6) =0 and 5a (a8) 7, 


The stationary point criterion. If (a,b) is an 
extremum of a well-behaved function f(z, y), and is not 
on the boundary of the domain of f(z, y), then (a,b) isa 
stationary point of f(x,y). (The converse statement is 
not necessarily true. ) 


(a, 6) 


4. A saddle point is a stationary point which is not an 
extremum. 


5. The AC — B’ criterion. Suppose (a,b) is a 
stationary point of f(x,y), and let 


arf 


maximum if A < 0, 
minimum if A > 0; 


(i) if AC — B? >0, (a,b) isa { 


(ii) if AC — B? <0, (a,)) is a saddle point; 


(iii) if AC — B? =0, the criterion gives no information. 


Unit 26 Vector calculus 


Section 1: Scalar and vector fields 


1. A scalar field function, or scalar field, is a 
function which to each point in a given region of space 
associates a unique scalar. The region is the domain of 
the scalar field. 


2. The contour curves (or contours) of a 
two-dimensional scalar field ¢(z,y) are the curves 
¢(x,y) = constant. The contour surfaces of a 
three-dimensional scalar field ¢(z, y, z) are the surfaces 
(x,y, z) = constant. 


3. A vector field function, or vector field, is a 
function which to every point in a given region of space 
(or of the plane) associates a unique vector. 


4. The field lines of a vector field are the family of 
continuous curves such that at every point the tangent to 
the curve is along the direction of the vector field at that 
point. 


Section 2: Differentiating scalar fields 


1. The symbol V, defined by 
O O O 
V =i—+j-—+k— 
"Ox oa, = Oz’ 
is a differential vector operator. 


2. For a three-dimensional scalar field $(z, y, z), the 
vector field 


Od. _ 9d, | OG 
— — —k 
Ox" 2 By? = Oz 

is called the gradient of ¢ and is denoted by grad ¢ or 


Vo. 


3. The direction of grad ¢ at a point P in the domain of 
the scalar field ¢ is perpendicular to the contour curve or 
surface passing through P. This direction is the normal 
to this contour at P. 


4. The vector field grad @ gives the magnitude and 
direction of the maximum (spatial) rate of change of the 
scalar field ¢ at any point P in the domain of the scalar 
field. 


5. The (spatial) rate of change of a scalar field ¢ at a 
point P in the direction of a unit vector e is given by 


e - (grad @). 


Section 3: The scalar line integral 


1. The scalar line integral of a vector field F along a 
curve C, called the path of integration, with 
end-points A and B, is defined as 


[Fae = tim SOF, or, 
Sc 


t= 4 
where the curve is divided into n segments with 
end-points P,(= A), Po, P3,..., Prii(= B); F; means 
os 


F(P;); 6r; means P;P;+1; and the limit is taken in such a 
way that n becomes very large and the segments all 
become very short. The end-points A and B may 
coincide, in which case the path of integration is closed; 
if they do not coincide, it is open. 


2. A three-dimensional curve can be expressed in terms 
of a parameter t by the parametric equations 


c=a2(t), y=y(t), z=2(t) (a<t<b). 


3. To evaluate a scalar line integral / F - dr: 
e 


(i) describe the path C by means of a set of parametric 
equations; 
(ii) use these parametric equations and the formula 
ae dz. ay ae 
eee ee 
to express the vector field F and dr/dt in terms of 
the parameter t, and hence rewrite the line integral 
as an integral with respect to ft: 


: dr 
[e-a= [ PO. Sat 
- . dt 


(iii) evaluate this integral with respect to t. 


4. The line integral of grad ¢ is independent of the path 
apart from its end-points: 


(grad ¢) « dr = 4(B) — (A). 


AB 


5. A vector field F is conservative if for every closed 
path C’ in the domain of F we have 


 F-dr=0 
C 


where the symbol means a line integral round the 


'g, 
closed path C’.. Otherwise the vector field is called 
non-conservative. 


Section 4: The curl of a vector field 


1. For a vector field F we define the component of 
curl F at a point P in the direction of a unit vector n to 
be 

1 


n- (curl F) = iim if F - dr, 
= C 


where C’ is a small closed curve in the plane through P 
perpendicular to n, which contains P and encloses an 
area A. (The direction of integration round C is chosen so 
that a screw turned in this direction advances along n.) 


2. The curl of a vector field F(z, y, z) = Fui+ Foj+ Fsk 
can be evaluated by using the derivative formula 
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an VS POrs 20 Fs) ; OF, OF3\. 
cul =VxP= (F i+ (F ae) 4 


3. The curl of a vector field F is a detector for rotational 
motion; for instance, the curl of a fluid velocity field 
describes the spin of a small float placed in the fluid. 


4. If F is a conservative vector field: 
(i) f F - dr = 0 for every closed path C in the domain 
c 
of F; 


F - dr is independent of the 


AB 
path taken between any two points A and B; 


(iii) curl F = 0; 
(iv) there exists a scalar field ¢ such that F = grad @; 


(ii) the line integral 


(v) / F - dr = $(B) — ¢(A), where F = grad ¢. 
AB 


5. Provided that the domain of F is sufficiently simple in 
shape, the converse of property 4(iii) above is also true, 
i.e. the equation curl F = 0 implies that F is conservative. 


Unit 27 Multiple integrals 


Section 1: The surface integral 


1. Suppose that f is a two-dimensional scalar field whose 
domain includes a region S of the (x, y)-plane. We divide 
the region S into a number N of area elements so that 
the zth element has area 6A;, and denote some point in 
the ith element by P;. Then the surface integral of f 
over the region S is defined to be 


N 
/ fad =lim)  f(P)6Ai, 
S 


7—4 
where the limit is taken in such a way that N becomes 
very large and the sizes of all the elements become very 
small. We call S the region of integration. 


2. The following procedure can be used to evaluate 
surface integrals. 


(i) Draw a diagram showing the region of integration, 
including the equations of its boundary lines. 


(ii) Mark on the diagram the minimum and maximum 
values of x, which we call here a and b respectively. 


(iii) Draw a typical strip of the region parallel to the 
y-axis and denote the formulas for the minimum 
and maximum values of y (which depend on the 
x-coordinate of the strip) by a(x) and G(r) 
respectively. 


(iv) The surface integral can now be written in terms of 
two single integrals: 


z=b y=A(ax) 
[sas / Flew) ay dz. 
S Ft yra(z) 


46 


(v) Evaluate the ‘inner’ integral, holding x constant: 


y=A(z) 
g(x) = 7 f(x,y) dy. 


=a(z) 


(vi) Evaluate the ‘outer’ integral: 


gah 
/ fdA= / g(x) da. 
S eae 


An analogous procedure, using strips parallel to the 
z-axis, can also be used in which we first integrate 
over x and then integrate over y. 


Section 2: Multiple integrals useful in 
mechanics 


1. Suppose that f is a three-dimensional scalar field 
whose domain includes a region B of space. The volume 
integral of f over the region B is defined as 


N 
/ fdV =lim )~ f(P:)6Vi, 
B i=l 


where B is divided into a number N of volume elements, 
the ith of which has volume 6V; and contains the point 
P;, and the limit is taken in such a way that all the 
volume elements become very small. We call B the 
region of integration. 


2. If p is the scalar field giving the local density of a 
body occupying a region B, then the mass of the body is 


m= | pay. 
B 


The moment of inertia of this same body about some 
given axis is 


be [ (od) dV, 


where d denotes distance from this axis. 


3. The mass of a hollow sphere of uniform density p with 
external radius a and internal radius 6 is <rp(a° — 6), 
and its moment of inertia about an axis through its 
centre is = mp(a° — b°). The corresponding formulas for a 
solid sphere are obtained by setting b = 0. 


4. To evaluate a surface integral / f dA using polar 


Ss 
coordinates (r,0) the following procedure can be used. 


(i) Draw a diagram showing the region of integration, 
S, and its bounding lines. 


(ii) Write the surface integral in terms of r and @ using 
the formulas 


x= rT cos8, 
y=rsin8, 
6A =r 6r 60. 


In the limit of very small area elements, this last 
formula becomes dA = r dr dé. 


(iii) Find the minimum and maximum values of 6 in S, 
say 9, and 62. 


(iv) For a strip of fixed 0, indicate the minimum and 
maximum values of r, say a and b (here assumed 
independent of 6, for simplicity). 


(v) Write the surface integral in terms of two single 
integrals: 


6=82 H==0 
/ / f(r, 9)r | dé. 
6=61 r=4 


(vi) Evaluate the inner integral, then the outer integral. 


Section 3: The volume integral 
The following procedure can be used to evaluate a 


f dV using Cartesian coordinates. 

B 

(i) Draw two diagrams, showing (a) the region B of 
integration and (b) the projection S of this region on 
the (x, y)-plane. (The projection of a point (z, y, z) 
on the (z, y)-plane is the point (x,y), and the 
projection of B comprises the projections of all the 
points in B.) 


volume integral 


(ii) Within the region B, draw a column parallel to the 
z-axis and mark its end-points, which are the 
equations of the top and bottom surfaces of B 
written in the form z = y(z,y) and z = d(z,y) 
respectively. 


(iii) Evaluate the single integral 


z=y(az,y) 
sag = i f(x,y, z) dz. 


=$(zx,y) 
(iv) Evaluate the surface integral of g(x,y) over the plane 


region S. The result is : f dV. 
B 


Similar procedures, with different orders of integration, 
can also be used. 


Section 4: Changing variables in three 
dimensions 


1. The cylindrical polar coordinates (r,0,z) of a 
point P in space are related to its Cartesian coordinates 


(x,y,z) by 


xz =rcos8, 
y=rsin8, 
22, 


where r > 0, 0 < 0 < 2m, —c0 < z < oo. That is, (r, 6) are 
the plane polar coordinates of the projection of P onto 
the (z, y)-plane. 


The volume element in cylindrical polar coordinates is 
6V =r 6r 60 6z. 
This becomes dV = rdrd@dz in the limit. 


2. To evaluate a volume integral using cylindrical polar 
coordinates, a procedure similar to the one given above 
for Cartesian coordinates can be used; the only difference 


is that we use plane polar coordinates r,@ in place of 
plane Cartesian coordinates x,y throughout. 


3. The spherical polar coordinates (r,0,@) of a point 
P are related to its Cartesian coordinates (x,y, z) by 


x =rsin@cos @¢, 

y =rsin@sin ¢, 

£208 0; 
where r > 0,0< ¢ < 27,0 <6< 7. In this case r is the 
distance from P to the origin O, @ is the angle between 
the line OP and the positive z-axis, and ¢ is the angle 


between OQ and the positive z-axis, where @ is the 
projection of P onto the (z, y)-plane. 


The volume element in spherical polar coordinates is 
6V =r’ sin 0 ér 50 4¢. 
This becomes dV = r? sin 6dr d6 d¢ in the limit. 


4. To evaluate a volume integral using spherical polar 
coordinates, use a procedure similar to the one given for 
Cartesian coordinates by expressing the volume integral 
as three repeated single integrals. 


Unit 28 Moments and circular motion 


Section 1: The components of a force 


1. The component of a force F in a direction making an 
angle a with that of the force is equal to F' cosa, where 
= 


2. If a particle is static (in equilibrium), then the sum of 
the components of the forces acting on the particle in any 
direction is equal to zero. 


Section 2: The moment of a force 


1. The magnitude [I of the moment of a force about a 
fixed point O is the product of the magnitude F' of the 
force and the perpendicular distance d from the point O 
to the line of action of the force, that is, 


P= Fe 


2. The total moment of the forces acting in two. 
dimensions on a body about any point is the sum of the 
anticlockwise moments minus the sum of the clockwise 
moments. 


3. If a rigid body is static (in equilibrium), then both the 
total force on the body and the total moment of forces 
acting on the body about any point are zero. 


AT 


Section 3: Motion in a circle 


1. In terms of polar coordinates, the radial and 
transverse components of the position, velocity and 
acceleration of a particle moving in a circle, r = constant, 
are as follows. 


Radial Transverse 
component | component 


Position 


Velocity 


Acceleration 


2. For circular motion with radius r and angular 
velocity w = 6, the velocity of the particle is v = rw, 
tangential to the circle. Both v and w are positive for 
anticlockwise motion. The component of acceleration 
towards the centre of the circle is rw? = v?/r. For 
non-uniform motion, the acceleration also has a 
transverse component rw. 


3. For uniform circular motion, with radius r and 
constant angular velocity w, the acceleration has 
magnitude rw” = v’/r and is directed towards the centre 
of the circle. The time taken for one complete revolution 
of the circle is T = 277/|w}. 


Section 4: Using vector notation 


1. The torque I, about a fixed point O, of a force F 
acting at a point which has position vector r relative to 
the point O, is defined to be 


=x F. 


2. The torque of all the gravitational forces acting on a 
system of particles or extended body can be evaluated by 
considering the total weight of the system to be 
concentrated at its centre of mass. 


3. If asystem of particles or extended body is in 
equilibrium, then the total external force is equal to zero, 
as is the total torque of the external forces about any 
fixed point. 


4, With respect to plane polar coordinates |r, 6], the unit 
vectors e, and eg (in the radial and transverse directions, 
respectively) are given by 


e, = cos#i+sin@j, eg = —sin@i+cos6@j. 
The vectors e,,eg9,k form a right-handed triad of unit 


vectors. The directions of e, and eg vary with time, so 
that 


e,=0eg and eg = —@e,. 


5. The position, velocity and acceleration of a particle 
moving in a circle of constant radius r in the (2, y)-plane, 
with centre at the origin, are 


= Te. 
v=r= reg 


and a= f= 0 6. + roep. 
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6. The angular velocity w for a particle in circular 
motion is a vector whose magnitude is equal to the 
angular speed, and whose direction is along the axis of 
rotation, in the sense in which the rotation would drive a 
right-handed screw. If the unit vector k is perpendicular 
to the plane of motion and the centre of the circle lies on 
the z-axis, then 


w — 6k. 


7. For circular motion with angular velocity w, the 
velocity of a particle is given by 


P= xy; 
where r is the position vector of the particle relative to 
an origin on the axis of rotation. 


Unit 29 Angular momentum and rigid 
bodies 


Section 1: Angular momentum 


1. The strong form of Newton’s third law states 
that for a pair of interacting particles, the force Fi2 on 
Particle 1 due to Particle 2 is equal in magnitude but 
opposite in direction to the force F2; on Particle 2 due to 
Particle 1, that is, 


Fi2 = =F ai, 


and each of these forces acts along the line which joins 
the two particles. 


2. For a particle which has position vector r relative to 
an origin O and linear momentum p = mr, the angular 
momentum I relative to O is 


f=T «p= Tf <a 


3. The torque law for a single particle states that the 
rate of change of the particle’s angular momentum 
relative to a fixed point is equal to the applied torque 
relative to that point, that is, 


eS 
In the special case where the total torque acting on the 
particle is zero, the angular momentum of the particle 
about the fixed point is constant. In other words, if 7 = O 
then lI is constant. This is the law of conservation of 
angular momentum for a single particle. 


4. The total angular momentum L of a system of n 
particles, relative to a given origin, is the vector sum of 
the angular momenta of the individual particles, that is, 


Tm nr 
= » ki = Pe x EES, 
i=1 i=1 
where r; is the position vector of the ith particle and m; 
is its mass. 


5. The total external torque I acting on a system 
of n particles, relative to a given origin, is defined to be 


the vector sum of the external torques which act upon 
the individual particles, that is, 


re*' ome ext 
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6. The torque law for a system of particles states that 
the rate of change of the total angular momentum about 


_ the origin is equal to the total external torque about the 


origin, that is, 

L _ rex! 
This law applies to any system for which the inter-particle 
forces satisfy the strong form of Newton’s third law. 


In particular, if the total external torque is zero then the 
total angular momentum L is constant. This is the law 
of conservation of angular momentum for a system 
of particles. 


Section 2: Rigid-body rotation 


1. A rigid body is a many-particle system with the 
property that all the inter-particle distances remain 
constant in time. The inter-particle forces in a rigid body 
are assumed to satisfy the strong form of Newton’s third 
law, so that the torque law L=rt applies. 


2. Rigid-body rotation about a fixed axis is a 
rigid-body motion for which particles on the axis remain 
fixed, while other particles move in circles centred on and 
perpendicular to the axis. If the axis of rotation is chosen 
to be the z-axis, then all particles have the same angular 
velocity w = 6k. If r; is the position of Particle i, then its 
velocity is given by 


fo Ss ee TE 


3. The z-component of the total angular momentum L of 
the body described above is 


L. = 16, 


where 


i=] 


is the moment of inertia of the body about the axis of 
rotation (the sum being taken over all particles of the 
body). It follows from the torque law that 


Bee ee eee 
The kinetic energy of the body is given by 

T= 116. 
4. The moment of inertia of a continuous rigid body 
about a given axis of rotation is 


= / p(r)d?(r) aV, 


where p(r) is the density of mass at point r, d(r) is the 
distance of the point r from the axis of rotation, and the 
integration is over the volume of the body B. 


For a homogeneous body, the density p is constant. The 
moments of inertia of certain homogeneous rigid bodies 
about axes through their centres of mass are given in 
Table 1. 


Table 1: Moments of inertia of homogeneous rigid bodies 


Homogeneous rigid Diagram Dimensions Moment of inertia 
body of mass M; about given axis 
axis 


solid sphere; radius R 2 MR? 
axis through centre 
of mass 
inner radius a, 
outer radius R 
lan ' 


inner radius a, 
outer radius R 


hollow sphere; 
axis through centre 
of mass 


solid cylinder; 
axis through centre 
of mass and along 
axis of cylinder 


hollow cylinder; 
axis through centre 
of mass and along 

axis of cylinder 


radius R, 
length h 


solid cylinder; 
axis through centre 
of mass and 
perpendicular to 
axis of cylinder 
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Table 1 (continued) 


Homogeneous rigid Diagram Dimensions Moment of inertia 
body of mass M; about given axis 
axis 


hollow cylinder; inner radius a, | }M(R* +a’) + + Mh? 
axis through centre outer radius R, 

of mass and length h & + MR? + + Mh? 
perpendicular to ve ae 

axis of cylinder 


solid rectangular faces 
brick; perpendicular 


axis through centre to axis have 

of mass and sides of 
perpendicular to lengths a and b 
one pair of faces 


thin rod length h 
(cross-section of 

arbitrary shape); 

axis through centre 

of mass and 

perpendicular to 


rod 


5. The parallel axes theorem states that the moment where = 
of inertia of a rigid body about an axis A can be written .rel)2 
- . To =4 > milti™| 
as : 
<< 
2 
I=I¢+ Md, is the kinetic energy relative to the centre of mass. 
where M is the mass of the body, d is the distance of the (ii) The total angular momentum of the system about 


centre of mass G from the axis A, and Ig is the moment 
of inertia of the body about an axis B which passes 
through G and is parallel to A. 


the origin is given by 
L=Lc+Rx MR, 
a ——__where 
Section 3: Rotation about a moving axis with i= 3 rte! x mete! 
fixed orientation 3 
| is the angular momentum about the centre of mass. 
1. Consider a system of n particles, with centre of mass 
R, total mass M, total kinetic energy T and total angular 
momentum L about the origin, which is acted upon by a 


(iii) The total external torque of the system about the 
origin is given by 


t t 
total external force F°** corresponding to a total external Eo =i +R x,” 
torque I'°** about the origin. Suppose that the ith where 
particle has mass ms and position r;, and experiences an “ee . pel. ext 
external force F§*". Then the position of the particle = i ; 


relative to the centre of mass, eo is defined by ae 
rel is the total external torque about the centre of 
ry =ri —-R. mass. 
The centre of mass decomposition theorems are as 
follows. 2. For a rigid body rotating with angular velocity w = 6k 
(i) The total kinetic energy of the system is given by about an axis which passes through the centre of mass 


jee ee ; M|RI, and is always aligned with the z-axis, the kinetic energy 


and the z-component of the total angular momentum 
relative to the centre of mass are given respectively by 


To =110 and Lo-k=I6, 
where J is the moment of inertia of the body about the 


axis of rotation. 


3. For any system of particles, the rate of change of the 
total angular momentum relative to the centre of mass is 
equal to the total external torque relative to the centre of 
mass, that is, 


Le =Tc. 
This is the torque law relative to the centre of 
mass. 


4. If acylinder or sphere of radius R rolls (without 
slipping) along a rough plane, then the speed v of the 
centre is related to the angular speed |6| of the body 
about its centre by the rolling condition 


v = Rid}. 


Section 4: Rotation about an axis whose 
orientation changes 


1. If a body spins around an axis of symmetry (the spin 
axis), and this axis is fixed, then the total angular 
momentum of the body is given by 


L = J6n, 
where J is the moment of inertia with respect to the spin 


axis and @n is the angular velocity of the body about this 
axis. 


2. Ifa body spins around an axis of symmetry which is 
not fixed, and J,@n are defined with respect to the spin 
axis as above, then the approximation 


L~I/@n 


is valid provided that |f| is small compared with @. 


3. The angular momentum of a spinning top has a 
constant magnitude and a constant vertical component. 
It satisfies the equation 


L=akxL, 
where @ is the rate of precession of the spin axis about 
the vertical axis, and k is a unit vector in the vertically 
upward direction. If the top has mass M, moment of 
inertia J and angular velocity @n about the spin axis, 
then the rate of precession is given by 

am ee 

I6 

where c is the distance of the centre of mass from the 
point of contact with the ground. 


Unit 30 Planetary orbits 


Section 1: Review of energy and planar 
motion 


1. For particle motion in three or fewer dimensions, a 
force F acting on the particle is conservative if there 
exists a scalar field function U = U(r) = U(z, y, z) such 
that 


F = — grad U. 


ol 


The function U is the potential energy of the particle. 
In such a case, the total mechanical energy of the particle, 
E=imr’+U 


2 


(where r? = r-r = |r|”), is constant throughout the 


motion. 


2. In one dimension, the condition for a conservative 
force and the expression for the total mechanical energy 
become respectively 
dU 
F=—— and F= Img? + U, 
dx 


where F' is the z-component of the force and U = U(z). 


3. The unit vectors e,,eg corresponding to the plane 

polar coordinates [r, 6] (where x = rcos@, y = rsin@) are 

related to the Cartesian unit vectors i,j by the equations 
e, =cos#i+sin@j, eg = —sinfi+coséj 


and i=cos#ée,—sinfdeg, j=sinde, + cosOeg. 


4. The derivatives e,,e@9 are given in terms of e,,e9 by 


e,=0eg and eé9 = —@e,-. 


5. The position r, velocity r and acceleration fr of a 
particle moving in the (z, y)-plane are given in terms of 
plane polar coordinates by the equations 

r= 'e;, 

r=—re,+ ree, 
ld 


== reo. 


2 
ee —_— reas fg ‘: me 
r= (7—T0@ )e, + ae 


Section 2: Kepler’s laws and Newton’s law of 
gravitation 


1. Kepler’s laws of planetary motion are as follows. 


(i) Each planet moves in an ellipse, with the Sun at 
one focus. 


(ii) The line joining a planet to the Sun sweeps out 
equal areas in equal times. 


(iii) The square of the orbital period of a planet is 
proportional to the cube of the semi-major axis of 
its orbit. 


2. (i) The ellipse shown in the figure has semi-major 
azis a and semi-minor azis b. With respect to the 
axes O’ XY, its equation is 
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(ii) The area of the ellipse is 7ab. 


(iii) The focz of the ellipse are located a distance q to 
either side of the centre O’, where 


g= Va? =- 6. 


(iv) With respect to axes Oxy with origin at the 
right-hand focus, and in polar coordinates, the 
equation of the ellipse is 


l 
—=1+ecos8, 
¢ 
where 
b2 
[= — 
a 
is the semt-latus rectum, and 
q b? 
C223 = 4 
a a? 


is the eccentricity, where 0 <e < 1. 


3. Newton’s universal law of gravitation states that 
the gravitational force on a particle with mass m; and 
position vector r; due to another particle with mass m; 
and position vector r; is 


Fi; = ——_4;; — v1), 


where G is the gravitational constant. 


4. A body has a spherically symmetric mass 
distribution if (with the origin of coordinates at the 
centre of the body) its mass density function p(r) 
depends only on the radial distance 

lr] =r = f/x? + y? + 2?. The gravitational force exerted 
by such a body B on an exterior particle with mass m 
and position vector r is given by 


GmM 
F(r) = — grad U(r) = — 2 ers 
where e, = r/r is the unit vector in the outward radial 
direction, 
GmM 
oe =- 


is the gravitational potential energy of the particle, and 


M = | p(r") av 


is the total mass of the body. 


Section 3: Orbits 


1. A central force F, acting from the origin on a 
particle with position vector r, has the form 


F(r) = g(r)r, 
where g(r) = g(z, y, z) is a scalar function of position. If 
the particle has mass m and is not acted upon by any 
other force, then it has constant angular momentum 
L =r x mr. If the unit vector k is chosen in the 
direction of L, then 


(i) the motion of the particle is confined to the 
(x, y)-plane; 

(ii) L=JL,k, where L, = mr7@ is constant. The 
gravitational force due to a spherically symmetric 
body whose centre is fixed at the origin is a central 
force. In this case, the constancy of L, is equivalent 
to Kepler’s second law. 


2. An isotropic central force has the form 
F(r) = f(r)er, 
where r = ,/z?2 + y? + 2? is the distance from the origin, 


and e, =r/r. For a particle of mass m acted upon by 
such a force, the radial equation of motion is 
E Le 

me = f(r) ax pe 
which may be integrated (analytically or numerically) to 
give r(t). The gravitational force is isotropic, with 
f(r) = —GmM/r?’, where M is the central mass and G is 
the gravitational constant. (It is assumed that M is much 
larger than m, in order that the central mass should 
remain fixed at the origin.) 


3. For the particle moving under gravity as described 
above, the total mechanical energy FE is constant. It may 
be expressed as 
(i) E = 1m(7? +176") + U(r), 
where U(r) = —GmM_/r is the (actual) potential 
energy of the particle; 


(ii) E = imi? +U"(r), 
where 
eff ; oe GmM 
Veen 2mr2 r 


is the effective potential energy. Orbits are 
unbound if E > 0, and bound if US! < E <0, 
where U4 is the minimum value of U ety) 

4. The equation of gravitational orbit for the particle 


described above is 


l 
—-=1+ecos8, 
r 
where the semi-latus rectum | is given by 
L? 
~ m2MG’ 


If 0 < e < 1 then the orbit is an ellipse with eccentricity 
e, showing that Kepler’s first law is satisfied. If e > 1 
then the motion is unbound. 
In terms of the eccentricity and the semi-latus rectum, 
the total mechanical energy is 
GmM 
fer). 


1 
2l 


5. Kepler’s third law is also a consequence of this 
Newtonian model, since the period T of an elliptical orbit 


and its semi-major axis a are related by the equation 


Ar? 


T*? =ka°, where k = MG’ 


Section 4: Central forces and conservation 
laws 


1. If a scalar function ¢(r) depends upon position only 
through the radial distance |r| = r = (x? + y* + z7)1/2 | so 
that ¢(r) = ¢(r), then the gradient of ¢ can be written as 
grad @ = af 
dr 
where e, =r/r is the unit vector in the outward radial 


direction. 


2. Hence if F is an isotropic central force, that is, 


F(r) = f(r)er 


for some function f, then F is conservative, with a 
potential energy function U(r) such that 


ee 
The motion of a particle in this force field is planar, and 
has constant total mechanical energy 


E = 1m(# +176) + U(r). 


Unit 31 Fourier analysis 


Section 1: What Fourier analysis is used for 


1. A function f(t) is said to be periodic, with period 
7’, if 

F() = fle eae ee Ft 2, +3... =); 
that is, the graph of f(t) consists of repeating blocks of 
length T. We usually take T as the size of the smallest 
such block, but any multiple of this also satisfies the 
definition of a period. 


2. If the response of a mechanical system to an input 
Fo sinwt is Asin(wt + ¢), we call A/F the relative 
response and ¢ the phase shift. The frequency 
response of the system is described by giving A and ¢ as 
functions of w. For a linear system (i.e. one that can be 
modelled by a linear differential equation whose input is 
the right-hand side of the equation and whose output is 
the unknown function), the principle of superposition (see 
Unit 6) enables us to find the output produced by any 
input that is a linear combination of sinusoidal functions. 
That is, if C1,C2,... are constants and F\(t), Fo(t),... 
are sinusoidal functions, and if the input F,(t) acting 

Tr 


alone produces output 2z,(t), then the input = C,F;,(t) 
Tm 
al | 


produces the output := Carty: 


r=! 


3. The sum of a number of sinusoidal terms of different 
frequencies is a non-sinusoidal function and in general is 
also non-periodic; but if the frequencies of the sinusoids 
are all multiples of w, then their sum will be periodic and 
will have period 27/w. In this case we call w the 
fundamental angular frequency. 


Section 2: Fourier series for periodic functions 


1. If f(t) is a periodic function with period T, then the 
Fourier series for f(t) has the form 


F(t)=M+ S (An sin nwt + B, cos nut), 
aol 
where w = 27/T is the fundamental angular frequency. 
The constants in this series are given by 


1 ft 
m=z s(t) at, 


-T /2 


ee hk. 
Awe z | f(Hyeinnetd - (n= 15,2, 3,...)5 


=T72 

5 Tl? 
B,== F(t) cosnat dt (n= 1, 2,3,....). 

T J _p/2 


2. A function is said to be continuous if its graph is an 
unbroken curve. A function is piecewise continuous if 
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its graph consists of continuous pieces, with only a finite 
number of such pieces on any finite part of the domain 
and each piece having a well-defined end-point. 


3. Theorem 2.3: Ifa periodic function and its 
derivative are piecewise continuous, then the function is 
equal to its Fourier series for all values of t at which there 
is no gap in the graph. 


Section 3: Some outstanding points 


1. A function f(t) is said to be even if 
f(t) = f(-—t) (for all t). 


The graph of an even function is symmetric under 


reflection in the vertical axis. The Fourier series for an 
even function contains no sine terms. 


An odd function f(t) is one satisfying 
f(t) = —f(-—t) (for all t). 


The graph of an odd function is symmetric under a 


rotation about the origin through 180°. The Fourier series 
for an odd function consists of sine terms only. 


2. We can often obtain useful approximations to a given 
periodic function f(t) by truncating its Fourier series, i.e. 
using a limited number of terms. 


Section 4: Extending the scope 


1. The theory described in this unit is not restricted to 
periodic functions of the variable t; any variable, for 
example xz, may be used. 


2. A non-periodic function whose domain is a finite 
interval of the t-axis can be represented by a Fourier 
series; we do this by extending the definition of the 
function to the entire t-axis in such a way that the 
extended function is periodic, and then computing the 
Fourier series for the extended function. This series then 
represents the original function on its own restricted 
domain. 


If the Fourier series is required to have sine (or cosine) 
terms only, it is necessary to use a periodic extension 
which is odd (or even). 


Unit 32 Partial differential equations 
Section 1: Introduction and definitions 


1. A partial differential equation is an equation 
involving one or more of the partial derivatives of some 
unknown function of two or more variables, and possibly 
also the value of the function itself and the values of the 
independent variables. (In contrast, an ordinary 
differential equation is a differential equation where 
the unknown is a function of just one variable, and which 
therefore does not involve partial derivatives.) Solutions 
of partial differential equations may be required to satisfy 
boundary conditions on the boundary of the domain 
of the unknown function. If one of the independent 
variables represents time, a condition imposed at the 
initial time is usually called an initial condition rather 
than a boundary condition. 
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2. An equation of the form 


0°U 0°U 0°U 
a(z,t) >> + b( tha ae c(x, t) > 
OU OU 
+d(z,t) 2 + e(z, t) a + f(x,t)U = g(z,t), 


where U(z,t) is the unknown function and at least one of 
the given functions a,b,c is not the zero function, is 
called a linear second-order partial differential 
equation. If in addition g(z,t) is the zero function, we 
say the equation is homogeneous; if the functions 
a,b,...,f are all constants, we say the equation has 
constant coefficients. 


3. Theorem 1: If wu and v are two solutions of any 
homogeneous linear partial differential equation, and A 
and B are arbitrary constants, then Au+ Bv is also a 
solution. 


Section 2: Properties of coiled springs and 
elastic strings 


1. A spring is in a state of tension if it is extended, 
and in a state of compression if it is compressed. The 
modulus of elasticity k of a spring is defined by 
K 

c= I” 
where K is its stiffness (denoted by k in Unit 7) and lo is 
its natural length. If a spring, or string, is in a state of 
tension, then for every point P along its length we define 
the tension T at P to be the magnitude of the force 
exerted by the part of the spring to one side of P on the 
part to the other side of P. The value of T is always 
positive for a spring or string in a state of tension. 


2. For a spring stretched between two fixed points A, B, 
we define the reference state as the configuration which 
makes the tension the same all along the spring. For a 
general state of the spring, still confined to the line AB, 
the tension T at a point P on the spring is given by 


ds 
T = (Io +k) — —k 
( 0 v ee ? 
where 
(i) Tp is the tension in the reference state, 
(ii) k is the modulus of elasticity, 


(iii) a2 is the position of P in the reference state, 
measured with respect to an x-axis lying along AB, 


(iv) s(x) is the position of P in the general state. 


If the spring is moving, so that s depends on the time t 


ds , 
as well as on 2, then = in the above formula is replaced 
z 


Os 


by 3a | 


ae ae 


3. A spring whose mass is not negligible is called a 
heavy spring. Vibrations of a stretched heavy spring in 
which all motion is along the direction of the spring itself 
are called longitudinal vibrations; vibrations in which 
the motion is at right angles to this direction are called 
transverse vibrations. Longitudinal vibrations can be 
modelled by the equation 

OU. og ; 

a? — 2 Of” @) 


where 


GQ) U Sex, 
(ii) c* = (To + k)/m, 


(iii) m is the mass of the spring per unit length in the 
reference state. 


If the ends of the spring are fixed at the points 
corresponding to zc = 0 and z = L, where L is given, the 
boundary conditions are 

Ui0,t) =e 

U(L,t) =0 } for all t. (2) 
Equation (1) is called the wave equation and the 
constant c in it is called the wave speed. 


Section 3: The method of separation of 
variables 


1. The method of separation of variables, applied to 
the wave equation (1) with boundary conditions (2), is to 
look for solutions of the form 


U(zx,t) = X(x)T(t), (3) 
where X,T are functions (of one variable) to be found. 


Substituting into the wave equation (1) and dividing by 
X(x)T(t) gives 


>, Ga 1 °C ) 
oY ee (4) 

X (zx) CT 
where p is a constant. So the function X satisfies the 
ordinary differential equation 

Xe) = ale). (5) 
The boundary conditions (2) imply 

¥ (0) = XG) — 6. (6) 


For p > 0 the only function X satisfying (5) and (6) is 
the zero function, leading to the trivial solution U = 0 of 
the original partial differential equation. For p<0a 
non-trivial solution is possible: 

X(z) = Fsinha; 
where pp = —h? (h > 0) and p is chosen so that 
sinhL = 0, i.e. 

nT 


Pa 
t oh 


where n is a positive integer. With this choice of w, when 
we have solved the equation 

* he a ce? uT 
implied by (4), the resulting solution U(z,t) = X(x)T(t) 
takes the form 


U Aad) = sin — (A. cos 


where a= 152.3... 


+ B,, sin 


niamct ewes ) 


2. A more general solution of the partial differential 
equation is (by Theorem 1): 


Uiz 4) = S - sin — (A. cos a + B, sin ee ees 2 


DL 


(oat | 
If the solution is required to satisfy initial conditions of 
the form 


U(x, 0) = p(x) 


= jnitial position of point z on spring, (8a) 


OU 
Sy (210) = ¥(2) 


= initial velocity of point x on spring, (8b) 


then the coefficients A, and By, can be found from the 
Fourier sine series expansions for the initial position and 
velocity functions ¢(x) and w(x) defined in Equations 
(8a) and (8b). 


Section 4: Transverse vibrations of a stretched 
string 


1. Transverse vibrations of a heavy spring or string can 
also be modelled by the wave equation, though the value 
of the wave speed c is not the same as for longitudinal 
vibrations. 


2. A boundary condition with the property that if two 
functions U; and U2 satisfy it then so do all linear 
combinations c;U; + c2U2 is said to be homogeneous; 
otherwise it is non-homogeneous. 


iti) 


3. The wave equation (1) with the non-homogeneous set 
of boundary conditions 


U(O,4)-= 0, (9a) 

U(L,t) = Rsinut, (9b) 
with w a constant, can be used to model transverse 
vibrations of a stretched string, one end of which is fixed 
and the other constrained to move sinusoidally. 
For this problem, separation of variables gives the 
particular solution 

Rs; 
pill a A a sin wt. 
sin(wL/c) 

The general solution is 

U(z, t) = Uo(z, t) © V(z,t), 
where V(z,t) stands for the right-hand side of Equation 
(7). If initial conditions are given as well, the constants 


Ay, By, can be determined using Fourier series, as in 
Section 3. 


Uo(z, t) — 
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Q Index of definitions 


A reference in this index beginning with the letter ‘p’ is to a particular page of this 


Handbook. 


A reference consisting of only numbers refers to the unit-by-unit summaries in 
Section 8 of this Handbook. Normally, the reference consists of two or three numbers, 
of which the first describes the unit, the second the section, and the third (where 
applicable) the numbered item, for those sections which are split into items. 

(For example, 29.3.4 means Item 4 of Section 3 of the entry for Unit 29.) 


A number on its own refers to the introductory part of the entry. (For example, 6 
means the Introduction of the entry for Unit 6.) 


absolute error 

absolute stability of numerical method 

absolute zero (of temperature) 

absolutely ill-conditioned 

absolutely well-conditioned 

AC — B? criterion 

acceleration (one-dimensional) 

acceleration (three-dimensional) 

acceleration vector 

addition of matrices 

addition of vectors 

air resistance 

amplitude 

angular frequency 

angular momentum 

angular velocity 

Argand diagram 

argument of complex number 

associated homogeneous equation of 
differential equation 

associativity 

augmented matrix 

auxiliary equation of differential equation 

auxiliary equation of recurrence relation 

axioms of algebra 


back substitution 
backward recurrence 
binomial theorem 
boundary condition 


Cartesian component of vector 
Cartesian form of complex number 
Cartesian unit vectors 

catenary 

central force 

centre of mass 

centre of mass decomposition theorems 
chain rule for ordinary differentiation 
chain rule for partial derivatives 
characteristic equation of matrix 
circular motion 

closure 

codomain 

coefficient of kinetic friction 
coefficient of restitution 

coefficient of static friction 
coefficients of differential equation 
coefficients of recurrence relation 
collision 

column matrix 


bd 728 
29.1.2, 29.1.4 
28.3.2, 28.4.6 
5.2 
5.3 


p. 9 
6.3, 32.1.1 


14.2.1 

5.1 

14.1.4 

18.3.3 
30.3.1, 30.4.2 
i eye 

29.3.1 

p. 14 

25.2.5 

rae oe: 
28.3.2, 28.4.5 
Bet 

p. 10 

15.3.4 

17.3.5 

15.3.4 

6 

1.1.3 

17.3.3 

20.1.1 


column vector 

combined heat transfer coefficient 
commutativity 

complementary function 

completing the square (for a quadratic) 
complex conjugate 

complex number 

component of force 

component of vector 

composite form for integration method 
composite function 

composition of functions 

conditioning of Gaussian elimination 
conduction 

conservation of angular momentum 
conservation of linear momentum 
conservation of (mechanical) energy 


conservative field 

conservative force 

consistent numerical method 
constant acceleration (motion with) 
constant catch strategy 

constant effort strategy 

constant function 

constant velocity (motion with) 
constant-coefficient differential equation 
constant-coefficient recurrence relation 
constraint 

continuous function 

contour 

contour curve 

contour surface 

convection 

convective heat transfer coefficient 
convergent numerical method 
critical damping 

cross product 

cubic (polynomial) 

curl 

cylindrical polar coordinates 


damped harmonic oscillator 
damped oscillations 
damped vibrations 
damping 

damping constant 

damping ratio 

dashpot 

dashpot constant 

dashpot force 


9.1 

28.1.1 

14.2.1 

18.4.7 

poli 

p. ll 

9.3.2 

ibe Be | 
29.1.3, 29.1.6 
Liga 

1.34, 3044 
30.3.3 

26.3.5 

rs Be ee 

19.3.1 

4.1.6 

3.4 


datum of potential energy 
de Moivre’s Theorem 
decomposition of matrix 
decomposition theorems 
degeneracy of linear programming 
problem 
degenerate vibrating system 
degree (of a polynomial) 
degrees of freedom 
dependent variable 
derivative 
derivative of vector function 
derived function 
determinant 
diagonal matrix 
diagonal of matrix 
difference of two squares 
differential vector operator 
differentiating partially 
differentiation 
direct iteration 


direction field of differential equation 


direction of vector 
displacement ratio 
displacement vector 
distinct eigenvalues 
distributivity 
domain 

dot product 


eccentricity (of ellipse) 
effective potential energy 
eigenvalue 

eigenvector 

elastic collision 

element of matrix 

ellipse 

empirical model 

equality of matrices 
error bound 

error function 

essential row interchange 
Euler’s formula for e’” 
Euler’s integration formula 
Euler’s method 


Euler-trapezoidal method 
even function 

explicit numerical method 
exploitation of populations 


exponential form of complex number 


exponential function of complex 
number 

exponential model for population 

external force 

external torque 

extrapolation 

extremum 


factorial 

factorization 

feasible region 

field line 

first-order differential equation 
first-order partial derivative 
first-order recurrence relation 


30.2.2, 30.3.4 
30.3.3 

21.1.1 

rake B 

17.3.4 

20.1.1 

30.2.2 

a 

20.1.2 

18.1.6, 18.4.3 
18.1.5 

9.3.1 

5.4 

18.4.4 

2.2, 6.5, 

oe Bey ie 
19.2.1 

g3.1 

19.1.4 

3.4 

5.4 


5.4 
3.4 

i733 
29.1.5 

3.1, 18.4.1 
25.4.2 


force of friction 
force of gravity 


force (one-dimensional) 

force (three-dimensional) 

force vector 

forced oscillations 

forced vibrations 

forward recurrence 

Fourier series 

Fourier’s law for conduction 
Fourier’s Theorem 

frequency 

frequency response 

friction 

frictional force 

frictional reaction 

function of a function rule 
function of many variables 
function of one variable 
fundamental angular frequency 
fundamental theorem of algebra 
fundamental theorem of calculus 


Gaussian elimination 

general solution of differential equation 
general solution of recurrence relation 
global error 

grad 

gradient of scalar function 

graph of function 

graphical method 

gravitational constant 

gravitational force 


gravitational potential energy 
gravity (force of) 


heat transfer 

heat transfer coefficient 

homogeneous boundary condition 
homogeneous differential equation 
homogeneous linear algebraic equations 
homogeneous recurrence relation 
Hooke’s law 

hyperbolic cosine function 

hyperbolic sine function 


identity element 

identity matrix 

image of a function 

imaginary (complex) number 

imaginary part of complex number 

implicit differentiation 

implicit numerical method 

independent variable 

induced instability of Gaussian 
elimination 

inelastic collision 

inhomogeneous differential equation 

inhomogeneous recurrence relation 

initial condition for differential equation 

initial condition for recurrence relation 

input-output principle 

integrating factor 
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15.3.4 
4.3.2, 28.4.2, 
30.2.3 
4.2.2 
15.1.4 
15.1.4 
8.2.2 
2:1 
1.4.2 
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1Z.2.1 
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15.3.4 
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p. 14 
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5.5 
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9.1 
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26.2.2 

20.2.2 
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10.2 
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L.12 

3.1 
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integration 

integration by parts 
integration by substitution 
internal force 
inter-particle force 
interpolating polynomial 
interpolation 

interval of absolute stability 
inverse (element) 

inverse iteration 

inverse matrix 

isotropic force 


Kepler’s first law 

Kepler’s laws 

Kepler’s second law 

Kepler’s third law 

kinetic energy 

kinetic energy of rotating rigid body 
kinetic friction 


law of gravitation 

length of curve 

limit of a function 

line integral 

linear combination of functions 
linear combination of matrix rows 
linear damping 

linear differential equation 
linear function 

linear momentum 

linear polynomial 

linear programming problem 
linear recurrence relation 
linear system of differential equations 
linear transformation 

linearity rule for integration 
linearly dependent matrix rows 
linearly independent functions 
local maximum 

local minimum 

local truncation error 

logistic curve 

logistic equation 

logistic model 

longitudinal vibrations 

lower triangular matrix 

LR method for eigenvalues 

LU decomposition 
lumped-parameter model 


magnitude of vector 
main diagonal of matrix 
matrix 

matrix addition 

matrix element 

matrix equality 

matrix multiplication 
matrix of coefficients 
matrix subtraction 
maximum 

Mean Value Theorem 
mechanical energy 
minimum 

Mixed Derivative Theorem 
modelling (mathematical) 
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p. 19 
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18.4.2 
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19.3.5 
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21.2:3 

20.4.3 
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30.2.1, 30.3.4 
30.2.1 

30.2.1, 30.3.1 
30:21. 30.3.5 
7.3.4 

29.2.3 

15.3.4 


30.2.3 
18.3.2 
p. 1l 
26.3.1 
6.1 
952:2 
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modified inverse iteration method 
modulus of complex number 
modulus of elasticity 
modulus of vector 

moment of force 

moment of inertia 
momentum (angular) 
momentum (linear) 
mortgage problem 
multiplication of matrices 
multipliers 


natural length of spring 

nested multiplication for polynomials 

Newton—Cotes formula 

Newton—Raphson method 

Newton’s first law 

Newton’s law of restitution 

Newton’s second law (one-dimensional) 

Newton’s second law (three-dimensional) 

Newton’s third law 

Newton’s universal law of gravitation 

non-homogeneous differential equation 

non-homogeneous recurrence relation 

non-singular matrix 

normal form of system of differential 
equations 

normal mode 

normal mode angular frequency 

normal reaction 

normal vector 


objective function 

odd function 

one-step numerical method 
optimal vertex 

optimum 

orbit 

order of differential equation 
order of recurrence relation 


parallel axes theorem 

parametric equation 

partial derivative 

partial differential equation 

partial differentiation 

partial fractions 

partial pivoting 

particle 

particular solution of differential equation 
particular solution of recurrence relation 
path of integration 

perfect dashpot 

perfect spring 

perfect square (polynomial) 

period 


periodic function 

phase (angle) 

phase shift 

phasor 

piecewise continuous 

pivot for Gaussian elimination 
pivot row 

plane polar unit vectors 

polar coordinates 


21.2.4 
5.3 

39. 9:4 

i434 

28.2.1 

27.2.9 29.953 
20.1 2.99.14 
17.3.1 

1.4.1 

20.2.1 

9.1 


Pat 
18.1.2 
18.4.7 
18.2.2 
4.2.1 
17.3.5 
4.2.3 
15.2.1 
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20.4.3 
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15.3.2 
14.3.9 


10.1 
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19.3244 
10.2 
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29.2.5 
26.3.2 

25.1.4, 25.3.1 
32:41:29 


9.3.1 
28.4.4, 30.1.3 
pO 


polar form of complex number 
polar unit vectors 

polynomial 

position (one-dimensional) 
position (three-dimensional) 
position vector 

potential energy 

power 

precession 

predictor-corrector method 


principal term in truncation error 


principal value of angle 
principal value range of inverse 
trigonometric function 
principle of superposition 
product rule for differentiation 
projection of vector 
proportionate birth rate 
proportionate death rate 
proportionate growth rate of 
population 


quadratic (polynomial) 
quadratic Taylor approximation 
quotient rule for differentiation 


radial component 

radiation 

range of projectile 

rate of precession 

rational function 

reaction 

real (complex) number 

real part of complex number 
recurrence relation 

region of integration 

relative error 

relative response 

relatively ill-conditioned 
relatively well-conditioned 
remainder for approximation 
repeated eigenvalue 
resonance 

restitution 

right-hand rule 

right-hand screw rule 
right-handed coordinate system 
rigid body 

rocket equation 

rolling condition 

root (of a polynomial) 

root of complex number 
rotation of rigid body 

row matrix 

row vector 

rule (of a function) 


saddle point 
scalar 

scalar field 

scalar line integral 
scalar product 
scale factor 
scaling of vector 
Schaefer model 
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28.4.4, 30.1.3 
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17.5.1 
29.3.4 
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9.0 
Pee, 29.552 
20.11 
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p. 10 


25.4.4 
14.1.1 
20.1.1 
26.3.1 
14.3.1 
1.3.3 
14.1.3 
3.3 


second derivative 

second-order differential equation 

second-order partial derivative 

second-order recurrence relation 

section function 

semi-latus rectum 

semi-major axis 

semi-minor axis 

separation of variables for first-order 
differential equation 

separation of variables for partial 
differential equation 

SI units 

simple harmonic motion 

simplex method 

Simpson’s integration formula 

Simpson’s method 

simultaneous linear algebraic equations 

singular matrix 

sinusoidal function 

size of matrix 

specific energy capacity 

specific heat 

speed 

spherical polar coordinates 

spring force 

spring potential energy 

spurious solution 

Square matrix 

stability of Gaussian elimination 

stability of numerical method 

static friction 

static rigid body 

stationary point 

steady-state vibration 

step length 

step size 

stiffness of spring 

strong damping 

substitution rule for integration 

subtraction of matrices 

subtraction of vectors 

sum rule for differentiation 

summation notation 

superposition 

surface 

surface integral 

symmetric matrix 


tangent approximation 
tangent plane 
Taylor approximation 


Taylor polynomial 


Taylor series 

Taylor series method 

Taylor’s Theorem 

temperature gradient 

tension of spring 

terminal speed 

thermal conductivity 

third derivative 

thrust of spring 

torque law about centre of mass 
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torque law about fixed point 


torque (of force) 

trajectory of projectile 

transfer coefficient 

transient 

transpose of matrix 

transverse component 

transverse vibrations 

trapezoidal integration formula 

trapezoidal method 

triangle inequality 

triangular matrix 

trigonometric functions 

trigonometric identities 

triple scalar product 

trivial constraint 

trivial solution of homogeneous 
linear equations 

truncation error 

turning points of motion 

two-step numerical method 


U value 

undamped angular frequency 
uniform circular motion 
uniform motion 

unit matrix 

unit vector 

universal law of gravitation 
upper triangular matrix 


vector 

vector addition 

vector differential operator 
vector field 

vector function 

vector product 

vector subtraction 

velocity (one-dimensional) 
velocity (three-dimensional) 
velocity vector 

vibration absorber 

volume integral 


wave equation 
wave speed 
weak damping 


yield of exploitation of populations 
zero function 


zero matrix 
zero vector 
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10 Notation 


The following symbols have standard meanings throughout the course. 


1 square root of —1 (see 5.1) 

g magnitude of acceleration due to gravity 
(approximately 9.81 ms~?) (see 4.3.2) 

sin, cos, tan p. 8 

sec, cosec, cot p. 8 


arcsin, arccos, etc. p. 9 


f(x) dx definite integral of function f(x) fom =a tor=b 
= belongs to 
[a, DJ the set of all z such thata<a<b 


Ox 


(x, y)-plane 


grad, curl see Unit 26 

lim pe ti 

log. x the logarithm of x to the base e = 2.71828... 

~ approximately equals 

! factorial 

a/be a~(bxc) 

fx the positive square root of x 

gin nth root of number x 

= Ye" 

|x| modulus of real number z (= Vx?) 

|z| modulus of complex number z 

z complex conjugate of complex number z 

|v| or v modulus or length of vector v 

Vv unit vector along v 

yle), fle) see p. 11 

Fz derived function of function f (p. 13) 

a derived function of f’ (p. 13) 

f<™ nth derivative of f (p. 13) 

first derivative, i.e. f’(x) where y = f(z) 

Zz = i.e. f’(t) where x = f(t) 

z ad i.e. f’(t) where x = f(t) 
es 

y an? ue f'(x) where y = f(z) 

e ory” second derivative, i.e. zs (#4) 

a"y or y\™) nth derivative 


indefinite integral of function f(z) 


the set of all x such thata<a2<b 
summation (see p. 5) 
x-axis of coordinate system 


plane containing z- and y-axes 
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